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Abstract 

We develop the general quantum measurement theory of non-Abelian anyons through in- 
terference experiments. The paper starts with a terse introduction to the theory of anyon 
models, focusing on the basic formalism necessary to apply standard quantum measure- 
ment theory to such systems. This is then applied to give a detailed analysis of anyonic 
charge measurements using a Mach-Zehnder interferometer for arbitrary anyon models. We 
find that, as anyonic probes are sent through the legs of the interferometer, superpositions 
of the total anyonic charge located in the target region collapse when they are distinguish- 
able via monodromy with the probe anyons, which also determines the rate of collapse. We 
give estimates on the number of probes needed to obtain a desired confidence level for the 
measurement outcome distinguishing between charges, and explicitly work out a number 
of examples for some significant anyon models. We apply the same techniques to describe 
interferometry measurements in a double point-contact interferometer realized in fractional 
quantum Hall systems. To lowest order in tunneling, these results essentially match those 
from the Mach-Zehnder interferometer, but we also provide the corrections due to pro- 
cesses involving multiple tunnelings. Finally, we give explicit predictions describing state 
measurements for experiments in the Abelian hierarchy states, the non-Abelian Moore- 
Read state at 1/ = 5/2 and Read-Rezayi state at u = 12/5. 
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1 Introduction 



One of the striking differences between two and three spatial dimensions is man- 
ifested in the allowed exchange statistics of quantum particles. In three spatial di- 
mensions, the allowed particle types can be classified according to irreducible rep- 
resentations of the permutation group. There are only two such one-dimensional 
representations, the trivial and alternating representations, for which exchange of 
any two particles introduces a factor of 1 and —1, respectively, to the wavefunction, 
corresponding to bosonic and fermionic statistics. Multi-dimensional representa- 
tions of the permutation group give rise to what is known as "parastatistics" [IJ, 
however, it has been shown that parastatistics can be replaced by bosonic and 
fermionic statistics, if a hidden degree of freedom (a non-Abelian isospin group) is 
introduced |l2l. 

In contrast, for two spatial dimensions the particles are classified according to rep- 
resentations of the braid group [3 |. The corresponding types of particles have been 
dubbed "anyons" II4I5II , and their exchange statistics are more precisely referred 
to as "braiding statistics." Braiding statistics described by multi-dimensional irre- 
ducible representations of the braid group [6 J give rise to non-Abelian anyons Q 

The group representation theory used to characterize particles' braiding statistics, 
however, becomes progressively cumbersome if one attempts to describe a sys- 
tem with several distinct "species" of anyons, especially those corresponding to 
multi-dimensional representations. Furthermore, one would typically like to con- 
sider systems in which there are processes that do not conserve particle number, a 
notion unsupported by the group theoretic language. To circumvent these shortcom- 
ings for systems with two spatial dimensions, one may switch over to the quantum 
field theoretic-type formalism of anyon models, in which the topological and al- 
gebraic properties of the anyonic system are described by category theory, rather 
than group theory. The structures of anyon models originated from conformal field 
theory (CFT) |7I8| and Chem-Simons theory [9J. They were further developed in 
terms of algebraic quantum field theory [llOlllll . and made mathematically rigorous 
in the language of braided tensor categories II12I13I14II . 

Surprisingly, even in our three-dimensional universe, there are physical systems 
that are ejfectively two dimensional and have quasiparticles - point-like localized 
coherent state excitations that behave like particles - that appear to possess such 
exotic braiding statistics. In fact, some of these are even strongly believed (though, 
thus far, experimentally unconfirmed) to be non-Abelian anyons. Anyon models 
describe the topological behavior of quasiparticle excitations in two-dimensional, 
many-body systems with an energy gap that suppresses (non-topological) long- 
range interactions, and hence an anyon model is said to characterize a system's 

^ In this paper, the term "anyon" will be used in reference to both the Abelian and non- 
Abelian varieties. 
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"topological order." 

The fractional quantum Hall effect is the most prominent example of anyonic sys- 
tems, so we will briefly review some relevant facts on the subject. (For a general in- 
troduction into the quantum Hall effect we refer the reader to Refs. [I15I16I17I18II .') 
The quantum Hall effect is an anomalous Hall effect that occurs in two dimensional 
electron gases (2DEGs) subjected to strong transverse magnetic fields (~ 10 T) at 
very low temperatures (~ 10 mK). Under these conditions, the Hall resistance R^y 
develops plateaus while R^x develops a series of deep minima as a function of 
the applied magnetic field. These plateaus occur at values which are quantized to 
extreme precision in integer [19] or fractional [|20ll multiples of the fundamental 
conductance quantum /h. These multiples are the filling fractions, usually de- 
noted V = Ng/N^ where A^e is the number of electrons and is the number of 
fundamental flux quanta through the area occupied by the 2DEG at magnetic field 
corresponding to the center of a plateau. At the plateaus, the conductance tensor is 
off-diagonal, meaning a dissipationless transverse current flows in response to an 
applied electric field. In particular, the electric field generated by threading an ad- 
ditional localized flux quantum through the system expels a net charge of ve, thus 
creating a quasihole. Consequently, charge and flux are intimately coupled together 
in the quantum Hall effect. 

In the fractional quantum Hall (FQH) regime, electrons form an incompressible 
fluid state that supports localized excitations (quasiholes and quasiparticles) which, 
for the simplest cases, carry one magnetic flux quantum and, hence, fractional 
charge ue. This combination of fractional charge and unit flux implies that they are 
anyons, due to their mutual Aharonov-Bohm effect. The fractional charge of quasi- 
particles in the v = 1/3 Laughlin state was first measured in 1995 [21]. Recently, a 
series of experiments asserting verification of the fractional braiding statistics has 
been reported [22.23.24.2512612711 . 

In the bulk of a FQH sample, the long-distance interactions between quasiholes are 
purely topological and may be described by an anyon model. Boundary excitations 
and currents of the Hall liquid are described by a 1 + 1 dimensional conformal field 
theory f2S\ whose topological order is the same as that of the bulk, when there is no 
edge reconstruction [I29I30II . These boundary excitations provide one way of cou- 
pling measurement devices to the 2DEG. A further connection between the physics 
of the bulk and CFT can be established following the observation in [|3TI that the 
microscopic trial wavefunction describing the ground state of the incompressible 
FQH liquid can be constructed from conformal blocks (CFT correlators). 

For the purpose of this paper, we are particularly interested in the possibility of non- 
Abelian statistics existing at several observed plateaus in the second Landau level 
(2 < u < 4), in particular u = 5/2, 7/2, and 12/5 II32I33I34I35 1. Predictions of 
non-Abelian statistics in FQH states originated with the paired state of Moore and 
Read [,31,1 , and was generalized by Read and Rezayi to a series of clustered non- 
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Abelian states ^M- At least for z/ = 5/2 (the Moore-Read state) and u = 12/5 (the 
k = 3, M = 1 Read-Rezayi state), these wavefunctions were found to have very 
good overlap with the exact ground states obtained by numerical diagonalization of 
small systems [,37.38il . 

Detailed investigations of the braiding behavior of quasiholes of the Moore-Read 
state were carried out in Ref. [39], and of the u = 12/5 state, as well as the other 
states in the Read-Rezayi series in Ref. [40]. Owing to the special feature of the 
Moore-Read state as a weakly-paired state of a p + ip superconductor of com- 
posite fermions 11411 . alternative explicit calculations of the non- Abelian exchange 
statistics of quasiparticles were carried out in the language of unpaired, zero-energy 
Majorana modes associated with the vortex cores II42I43II . (Unfortunately, this lan- 
guage does not readily adapt to give a similar interpretation for the other states in 
the Read-Rezayi series.) 

In addition to the proposed fractional quantum Hall states that could host non- 
Abelian anyons B 1136144 1. there are a number of other more speculative proposals 
of systems that may be able to exhibit non-Abelian braiding statistics. These in- 
clude lattice models [ 45.46.1 . quantum loop gases [47 48 49.50.1 . string-net gases 
[[51|52|53I54[[ . Josephson junction arrays [55], p + ip superconductors [[56I57I58[ [. 
and rapidly rotating bose condensates [[59I60I61[| . Since non-Abelian anyons are 
representative of an entirely new and exotic phase of matter, their discovery would 
be of great importance, in and of itself. However, as additional motivation, non- 
Abelian anyons could also turn out to be an invaluable resource for quantum com- 
puting. 

The idea to use the non-local, multi-dimensional state space shared by non-Abelian 
anyons as a place to encode qubits was put forth by Kitaev [45], and further de- 
veloped in Refs. [62.63 64 65 . 66^671681 . The advantage of this scheme, known as 
"topological quantum computing," is that the non-local state space is impervious 
to local perturbations, so the qubit encoded there is "topologically" protected from 
errors. A model for topological qubits in the Moore-Read state was proposed in 
Ref. [69], however braiding operations alone in this state are not computation- 
ally universal, severely limiting its usefulness in this regard. Nevertheless, one 
may still hope to salvage the situation by supplementing braiding in the Moore- 
Read state with topology changing operations [[7017 1[[ or non-topologically pro- 
tected operations fl2i to produce universality. The greater hope, however, lies in 
the A; = 3 Read-Rezayi state, for which the non-Abelian braiding statistics are es- 
sentially described by the computationally universal "Fibonnaci" anyon model (see 
Section [3. 6. 2[) . Consequently, the efforts in "topological quantum compiling" (i.e. 
designing anyon braids that produce desired computational gates) for this anyon 
model [[73I74I75[[ may be applied directly. 

The primary focus of this paper is to address the measurement theory of anyonic 
charge. This provides a key element in detecting non-Abelian statistics and cor- 
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rectly identifying the topological order of a system. Furthermore, the ability to 
perform measurements of anyonic charge is a crucial component of topological 
quantum computing, in particular for the purposes of qubit initialization and read- 
out. Clearly, the most direct way of probing braiding statistics is through experi- 
ments that establish interference between different braiding operations. In this vein, 
we will consider interferometry experiments which probe braiding statistics via 
Aharonov-Bohm type interactions [1761 . where probe anyons exhibit quantum in- 
terference between homotopically distinct paths traveled around a target, producing 
measurement distributions that distinguish different anyonic charges in the target. 
This sort of experiment provides a quantum non-demolitional measurement [ITTll 
and is ideally suited for the qubit readout procedure in topological quantum com- 
puting. 

The paper is structured as follows: 

In Section 121 we provide an introduction to the theory of anyon models, giving all 
the essential background needed to understand the rest of the paper, and establish- 
ing the connection with standard concepts of quantum information theory. 

In Section[3l we analyze a Mach-Zehnder type interferometer for an arbitrary anyon 
model. We consider a target anyon allowed to be in a superposition of anyonic 
states, and describe its collapse behavior resulting from interferometry measure- 
ments by probe anyons. We find that probe anyons will collapse any superpositions 
of states they can distinguish by monodromy, as well as decohere anyonic charge 
entanglement that they can detect between the target and outside anyons. We show 
how these measurements may be used to determine the target's anyonic charge 
and/or help identify the topological order of a system. We conclude this section by 
applying the results to a few particularly relevant examples. 

In Section m we consider a double point-contact interferometer designed for frac- 
tional quantum Hall systems. We give the evolution operator to all orders in tunnel- 
ing, and apply the methods and results of Section [3] to describe how superpositions 
in the target anyon state collapse as a result of interferometry measurements, and 
how to determine the anyonic charge of the target. We give detailed predictions 
for the Abelian hierarchy states, the Moore-Read state (u = 5/2,7/2), and the 
k = 3,M = 1 Read-Rezayi state (z/ = 12/5). 



2 Anyon Models 

In this section, we briefly review aspects of the theory of anyon models which are 
relevant to the rest of the paper. We follow the relatively concrete approach found 
in Refs. [ 4617811 . and develop some concepts in this formalism that are essential for 
the treatment of the measurement problem, such as the density matrix description 
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of states and the partial trace and partial quantum trace. 
2. 1 Fusion and Quantum Dimensions 

An anyon model has a finite set C of superselection sector labels called topologi- 
cal or anyonic charges. These conserved charges obey a commutative, associative 
fusion algebra 

axh = Y.K,c (2.1) 

cec 

where the fusion multiplicities N^^ are non-negative integers which indicate the 
number of different ways the charges a and h can be combined to produce the 
charge c. There is a unique trivial "vacuum" charge 1 G C for which N^-^ = Sac, 
and each charge a has a unique conjugate charge, or "antiparticle," a e C such that 

Nab = ^ba- (1 = i and i = a.) 

In order to have a non-Abelian representation of the braid group (details on braiding 
follow), there must be at least one pair of charges a and b in the theory which have 
multiple fusion channels, i.e. 

c 

The domain of a sum will henceforth be left implicit when it runs over all possible 
labels. Charges a which have J^c — 1 for every h must correspond to Abelian 
anyons (possibly bosons or fermions). Abelian and non-Abelian charges may also 
be distinguished by their quantum dimensions. The quantum dimension da of a 
charge a is a measure for the amount of entropy contributed to the system by the 
presence of a particle of type a. It may be found from the fusion multiplicities by 
considering the asymptotic scaling of the number of possible fusion channels for 
n anyons of charge a. For large n, this scales as d". Abelian charges have quan- 
tum dimension equal to 1, while non-Abelian charges quantum dimensions strictly 
larger than 1. The total quantum dimension of an anyon model is defined as 



1^ = Jj:dl (2.3) 



2.2 States, Operators and Inner Product 

To each fusion product, there is assigned a fusion vector space V^f, with dim V^, = 
N^^,, and a corresponding splitting space V^, which is the dual space. We pick 
some orthonormal set of basis vectors |a, b; c, pL) e V^^ ((a, b; c, //| e V^j^) for these 
spaces, where /i = 1, . . . , A^"^;,. If A^"^;, = 0, then V^'' is zero-dimensional and it has 
no basis elements. We will sometimes use the notation c E {a x b} to mean c such 
that A*"^^ 7^ 0. Allowed splitting and fusion spaces involving the vacuum charge 
have dimension one, and so we will leave their basis vector labels n — l implicit. 



6 



It is extremely useful to employ a diagrammatic formalism for anyon models. Each 
anyonic charge label is associated with an oriented line. It is useful in some contexts 
to think of these lines as the anyons' worldlines (we will consider time as increasing 
in the upward direction), however, such an interpretation is not necessary nor even 
always appropriate. Reversing the orientation of a line is equivalent to conjugating 
the charge labeling it, i.e. 



(2.4) 



The fusion and splitting states are assigned to trivalent vertices with the appropri- 
ately corresponding anyonic charges: 



{dc/dadbf^* 



aX\, b 



{a 



,b; c,(x\ 



abi 



(2.5) 



{djdjbf' "^YI " = \a, b; c, /i) G V:", (2.6) 




where the normalization factors (dc/dadb)^^^ are included so that diagrams are in 
the isotopy invariant convention throughout this paper. Isotopy invariance means 
that the value of a (labeled) diagram is not changed by continuous deformations, 
so long as open endpoints are held fixed and lines are not passed through each 
other or around open endpoints. Open endpoints should be thought of as ending on 
some boundary (e.g. a timeslice or an edge of the system) through which isotopy 
is not permitted. Building in isotopy invariance is a bit more complicated than just 
making this normalization change, but for the purposes of this paper, we can ignore 
the full details (which can be found in [I46I79II ). 

Inner products are formed diagrammatically by stacking vertices so the fusing/splitting 
lines connect 




I dadb 



d 



c 



(2.7) 



and this generalizes to more complicated diagrams. An important feature of this 
relation is that it diagrammatically encodes charge conservation, and, in particular, 
forbids tadpole diagrams. An important special case is c = 1, which shows that an 
unknotted loop carrying charge a evaluates to its quantum dimension 

= da = d-a. (2.8) 
The completeness relation for the identity operator on a pair of anyons with charges 
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a and b respectively is written diagrammatically as 




Any diagrammatic equation, such as this, is also valid as a local relation within 
larger, more complicated diagrams. Using Eq. (12.91) , Eq. (12.71 ) and isotopy, we get 
the following important relation, which expresses the compatibility of the quantum 
dimensions with fusion. 




(2.10) 



For general operators, we introduce the notation 



A, 



Ar. 



X 



V 



A' 



"T 

a; 



0,1 ,...,ar 



(2.11) 



where a capitalized anyonic charge label means a (direct) sum over all possible 
charges, so that the operator X is defined for acting on any n anyon input and 
m anyon output. The box stands for a linear combination of diagrams describing 
the action of the operator. The indices on operators will be left implicit when they 
are contextually clear (and unnecessary). Conjugation of a diagram or operator is 
carried out by simultaneously reflecting the diagram across the horizontal plane 
and reversing the orientation of arrows. 

Tensoring together two operators (on separate sets of anyons) is simply executed 
by juxtaposition of their diagrams: 



X 



Y 



(2.12) 
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2.3 Associativity 



The splitting of three anyons with charges a, 6, c from the charge d corresponds to 
a space V^'"' which can be decomposed into tensor products of two anyon splitting 
spaces by matching the intermediate charge. This can be done in two isomorphic 
ways 

yabc ^ yab ^ yec ^ yaf ^ ybc_ (^2.13) 

/ 

To incorporate the notion of associativity at the level of splitting spaces, we need 
to specify a set of unitary isomorphisms between different decompositions that are 
to be considered simply a change of basis. These isomorphisms (called F-moves) 
are written diagrammatically as 



ahc abc 




(2.14) 



d 



The same notion of associativity is, of course, true for fusion of three anyons. The 
associativity for fusion is given by F\ and together with unitarity, we have 



abcy 



j^abc 



(e,Q,/3)(/,M,i.) 



Fd ) 



-1 



(/,M,i/)(e,a,/3) 



(2.15) 



For fusion and splitting of more anyons, one does the obvious iteration of such 
decompositions. For this to be consistent, the F-symbols must satisfy a constraint 
called the Pentagon equation. One also imposes the physical requirement that fu- 
sion and splitting with the vacuum charge does not change the state. This means in 
particular that in diagrams, we may move, add, and delete vacuum lines at will (in 
fact, we already did this in Eq. 12.81) . Despite the constraints, the F-symbols have a 
certain amount of "gauge freedom," which comes from the fact that we are free to 
choose bases for the vertex spaces. 

We will also need the F-move with one of its legs bent down 



e 



/3 

^d 




(2.16) 



which is also a unitary transformation. From Eq. (12.91) . we immediately find that 

{Kb) 



■pab 
ab 



l(c,fj,,i/) 



1 1 'J/J.Ul 

O'aO'b 



(2.17) 



and more generally, applying Eqs. (12. 9|) and (12.71) gives a relation between the two 
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types of F-symbols 



:pab 
cd 



(e,a,/3)(/,/x,i/) 



dadd 



F 



ceb 



{a,a,ij.){d,l3,v) 



(2.18) 



2.4 Trace and Partial Trace 



The trace over operators formed from bras and kets is defined in the usual way. To 

translate the trace into the diagrammatic formalism, one defines the quantum trace, 
denoted Tr, by closing the diagram with loops that match the outgoing lines with 
the respective incoming lines at the same position 



TvX = Tr 



A-i A, 



X 



A' A' 



Alf. . . An 



X 



(2.19) 



Connecting the endpoints of two lines labeled by different anyonic charges violates 
charge conservation, so such diagrams evaluate to zero. The operator X e 
may be written as 



X, e K^- 



(2.20) 



(note that this decomposition is basis independent), which may be used to relate the 
quantum trace and the standard trace of bras and kets via 



TrX 



E|TrX„ 

c 



TrX = E dj:rX^ 



(2.21) 



Note that these are the same when the overall charge of the system is the vacuum 
charge (or any Abelian charge for that matter). 



We also need to define the partial traces for anyons. Since we have not yet intro- 
duced braiding, in order to take the partial trace over a single anyon B, the planar 
structure requires that it must be one of the two outer anyons (i.e. the first or last in 
the lineup). Physically, this corresponds to the fact that one cannot treat the subsys- 
tem excluding B as independent of B if this anyon is still located in the midst of the 
remaining anyons. The partial quantum trace over S of an operator X e yA'^'"'A'^B' 
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is defined by looping only the line for anyon B back on itself 



X 



"5 



A' A' 



(2.22) 



and for X E V^/'^-;'^? as 



Al An 



TtbX = B^' 



X 



A' A' 



(2.23) 



To relate the partial quantum trace to the partial trace, we implement factors for 
the quantum dimensions of the overall charges of the operator before and after the 
partial trace 



Tr^X, 



dr 



-fdf 



(2.24) 



where 



TrsX, = J2 [TrsX, 

/ 



TrsX, 



A' 



(2.25) 



The partial trace and partial quantum trace over the subsystem of anyons 

B = {Bi, . . . , Bn} that are sequential outer lines (on either, possibly alternating, 

sides) of an operator is defined by iterating the partial quantum trace on the B 

anyons 



Tr^ = TrR, . . . Tr 



Tre = Tr 



Tr 



(2.26) 



Iterating these over all the anyons of a system returns the trace and quantum trace, 
respectively, as they should. 



Using Eq. (12.171) and the fact that tadpole diagrams evaluate to zero, we may cal- 
culate the partial quantum trace for basis elements of two-particle operators 
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From this, we see that the partial trace acting on ket-bra elements seems to behave 
as the usual partial trace, and one might think it should be treated as such. However, 
things are a bit more subtle than this, since these bras and kets do not have the usual 
tensor product structure. When considering tensor products of operators, it is the 
partial quantum trace that behaves in the appropriate manner for a partial traces 
(i.e. as in the usual basis independent definition of partial trace), indicating that this 
should be treated as the usual notion of partial trace. Specifically, tracing over the 
set of anyons B on which the operator Y acts, we have 

Trs [X®Y] = XTrY (2.28) 

TVB [X y] = ^ Kb^aTrY,. (2.29) 



2.5 Braiding 



The unitary braiding operations of pairs of anyons, also called i?-moves, are written 
as 

Rab = , Rl = R'b' = V , (2.30) 

which are defined through their application to basis vectors: 



Rab\a,b;c,i^)^Y.\Rf] \b,a;c,u) (2.31) 

% = E K - 'Y. (2.32) 

cT 



and similarly for ^, which, by unitarity, satisfy 

For braiding to be consistent with fusion, it must satisfy constraints (the Hexagon 
equations), which essentially impose the property that lines may be passed over or 
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under vertices respectively (i.e. braiding commutes with fusion), and which imply 
the usual Yang-Baxter relation for braids. 



With the ability to braid, one also gains the ability to trace out any anyon in a 
system, not just those situated at one of the two outer positions of a planar diagram. 

To do so, one simply uses a series of braiding operations to move the anyon to 
one of the outside positions. In general, if one applies a different series of braids 
to move the anyon to the outside position before tracing, it will give a different 
outcome. Consequently, the braiding path applied to an anyon before closing its 
charge line should be included as part of the definition of the partial (quantum) 
trace. Physically, this corresponds to specifying the path (with respect to the other 
anyons) by which the traced anyon is removed from the system in consideration. In 
this paper, we will diagrammatically indicate the removal path of the anyon being 
traced out whenever the issue arises. 



The braiding matrices satisfy the ribbon property 



A 



//A 



c 



(2.33) 



where Oa is a root of unity called the topological spin of a, defined by 

ea = e-a = d-^'TvRaa = E | = ^ (^^^^ " ^^.34) 

When applicable, this is related to the (ordinary angular momentum) spin or 
CFT conformal scaling dimension of a, by 

= e*2-^^ (2.35) 

The topological -S-matrix is defined by 

Sab = V-^Tr [R,aRa,] = V-^ ^ ^-^-^^^^ = 3 aQO^ ' ^^'^^^ 

One can see from this that Sab = S^a — S^^, and da — Sia/Sn. A useful property 
for removing loops from lines is 



CD 

I 



(2.37) 



An anyon model is "modular" and corresponds to a TQFT (topological quantum 

field theory), if its monodromy is non-degenerate, i.e. for each a ^ 1, there is some 
b such that RfjaRab 7^ Ia6, which is the case iff the topological .S-matrix is unitary. 
For such theories, the S'-matrix, together with Tab = ^a^ab represent the generators 
of the modular group PSL (2, Z). 
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The monodromy scalar component 



M, 



ab 



Tr [RbaR. 



ubl 



TrI, 



ab 



dadb a 




SabSu 



SiaS- 



(2.38) 
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is an important quantity, typically arising in interference terms, such as those oc- 
curring in experiments that probe anyonic charge. It may also be written the 1,1 
component of the operator B'^ — F'^R^F. Since B"^ is a unitary operator, we must 



have|Mab| < 1. Indeed unitarity implies that when |Mab| 
of B'^ is non-zero, hence 




M, 



ab 



b 



1, only the 1, 1 element 



(2.39) 



so that the braiding of a and b is Abelian. The monodromy of a and b is trivial if 
Mab = 1. If ^ and |Mbe| = 1 for some e, then the relation 



(2.40) 



follows from the diagrammatic equation 




be e 




(2.41) 



2. 6 States and Density Matrices 



To describe the state of anyons in a system using a state vector, one must spec- 
ify all the splitting channels starting from vacuum. For example, in order to have 
anyons with charges a and b with overall charge c, one must also have a c charge 
somewhere, and one would write a general state of this form as 



'^a,b,c,^l\aMc,^J)\c,c^A) 

a,b,c,fi 

a b c 

a,b,c,iJ. {dadbdc) 

Unfortunately, describing states in this manner can become cumbersome, and, as 
usual, does not naturally accommodate the restriction to subsystems, so it is better 
for us to use the density matrix formalism. 

The density matrix for an arbitrary two anyon system is 




(2.42) 
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a,a',b,b' 

C,fJ,,fJ,' 



P{a,b,c,^i){a',b',c,fj.') ^JM 



a,a',b,b' {dadbda'dh'dl) 



1/4 ■ 



(2.43) 



The overall charge c must match up between the bra and the ket because of charge 
conservation. The normalization is chosen such that the trace condition for density 
matrices takes the form 

Tr [p] = P(a,b,c,M)(a,fe,c,/.) = 1 (2.44) 

a,b,c,fj, 

The factor 1/dc could, of course, be absorbed into P{a,b,c,fi)[a' ,b' ,cy) (as a matter of 
convention), but then the would appear in the summand of Eq. (12.441) . Further- 
more, the density matrix is written this way so that one can naturally think of it as 
p = Tr^ [p'] , the partial quantum trace over C of a density matrix that describes the 
actual entire system 



P = P(a,b,c,^i){a',b',c',fi') \0',b;c,p) \c,c; 1) (c, c'; 1| {a,b';c,p\ 



a,b,c,fM 
a' ,b' ,c' ,fj.' 



a b c 




C 



P{a.,b,c,fJ.){a',b',c',fi') ^2 

a,b,c,fi {dadbdcda'dh'dc')^^'^ 
a',b',c',ii' 




a' 



b' 



which only has vacuum overall charge. In other words, the entire system really 
has trivial total anyonic charge, but by restricting our attention to some subset of 
anyons, we have a reduced subsystem with overall charge c. Tracing over the C 
anyon (which imposes c = c') physically represents the fact that it is no longer 
included in the system of interest, and cannot be brought back to interact with the 
A and B anyons. Because of this, we are restricted to a subsystem which may only 
have incoherent superpositions of different overall charges c (i.e. one must keep 
track of the C anyon to allow access to coherent superpositions). The manifestation 
of this property in p is exhibited by the charge c matching in the bra and the ket (or 
diagrammatically as the charge c line connecting p and p'). The generalization to 
density matrices of arbitrary numbers of anyons should be clear. 

When considering the combination of two sets of anyons A = {Ai, . . . , Am} and 
B = {Bi, . . . , Bn}, the anyons of system A are unentangled with those of system 
B if the density matrix of the combined system is the tensor product (in some basis) 
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of density matrices of the two systems p^^ = ® . This essentially means the 
creation histories of the two different systems do not involve each other. There is a 
specific aspect of entanglement in anyonic systems that we will call anyonic charge 
entanglement, which is encoded in the anyonic charge lines connecting anyons. 
The systems of anyons A and B are said to have no anyonic charge entanglement 
between them if p^-^ G V^^'"''^r ®V^,^'"''^r . This is represented diagrammatically 
as being able to write the combined state such that there are no non-trivial charge 
lines connecting the anyons of A with those of B. 



3 Mach-Zehnder Interferometer 

In this section, we consider, in detail, a Mach-Zehnder type interferometer II80I81II 
(see Fig.[T]) for quasiparticles with non-Abelian anyonic braiding statistics, extend- 
ing the analysis begun in Ref. [|82ll . This will serve as a prototypic model of inter- 
ferometry experiments with anyons, and the methods used in its analysis readily 
apply to other classes of interferometers (e.g. the FQH double point-contact in- 
terferometer considered in Section H]). This interferometer was first considered for 
non-Abelian anyons in Ref. If83l . but only for anyon models described by a discrete 
gauge theory-type formalism in which individual particles are assumed to have in- 
ternal Hilbert spaces, and which use probe anyons that are all identical and have 
trivial self -braiding. Unfortunately, this excludes perhaps the most important class 
of anyon models - those describing the fractional quantum Hall states - so we 
must dispense with such restrictions. We will abstract to an idealized system that 
supports an arbitrary anyon model and also allows for a number of desired manipu- 
lations to be effected. Specifically, without concern for ways to physically actualize 
such manipulations, we posit the experimental abilities to: (1) produce, isolate, and 
position desired anyons, (2) provide anyons with some manner of propulsion to 
produce a beam of probe anyons, (3) construct lossless beam-splitters and mirrors, 
and (4) detect the presence of a probe anyon at the output legs of the interferometer. 

The target anyon A is the composite of all anyons Ai,A2,... that are located in- 
side the central interferometry region, and so may be in a superposition of states 
with different total anyonic charges. Since these anyons are treated collectively 
by the experiment, we ignore their individuality and consider them as a single 
anyon A capable of charge superposition. We will similarly allow the probe anyons, 
Bi, . . . , B]^ to be treated as capable of charge superposition (though this would cer- 
tainly be more difficult to physically realize). The probe anyons are sent as a beam 
into the interferometer through two possible input channels. They pass through a 
beam splitter Ti, are reflected by mirrors around the central target region, pass 
through a second beam splitter T2, and then are detected at one of the two possible 
output channels by the detectors Dg. When a probe anyon B passes through the 
bottom path of the interferometer, the state acquires the phase e*^', which results 
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D 



Fig. 1. A Mach-Zehnder interferometer for an anyonic system. The target anyon(s) A in the 
central region shares entanglement only with the anyon(s) C outside this region. A beam of 
probe anyons Bi, . . . , Bn is sent through the interferometer, where Tj are beam splitters, 
and detected at one of the two possible outputs by Ds- 



Fig. 2. The transmission and reflection coefficients for a beam splitter. 

from background Aharonov-Bohm interactions |[76l . path length differences, phase 
shifters, etc., and is also acted upon by the braiding operator Rba^ which is strictly 
due to the braiding statistics between the probe and target anyons. Similarly, when 
the probe passes through the top path of the interferometer, the state acquires the 
phase e*^" and is acted on by R~2b. 

Using the two-component vector notation 

to indicate the direction (horizontal or vertical) a probe anyon is traveling through 
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the interferometer at any point, the lossless beam splitters (S4] (see Fig. ^ are 
represented by 



T- 

^3 



(3.2) 



1 2 1 2 

(for J = 1, 2), where \tj\ + \rj\ = 1. We note that these matrices could be mul- 
tiplied by overall phases without affecting any of the results, since such phases are 
not distinguished by the two paths. 



When considering operations involving non-Abelian anyons, it is important to keep 
track of all other anyons with which there is non-trivial entanglement. Indeed, if 
these additional particles are not tracked or are physically inaccessible, one should 
trace them out of the system, forgoing the ability to use them to form coherent 
superpositions of anyonic charge. We assume that the target anyon has no initial 
entanglement with the probe anyons, so their systems will be combined as tensor 
products, with no non-trivial charge lines connecting them before they interact in 
the interferometer. 



The target system involves the target anyon A and the anyon C which is the only 
one entangled with A that is kept physically accessible. Recall that these anyons 
may really represent multiple quasiparticles that are being treated collectively, but 
as long as we are not interested in operations involving the individual quasiparticles, 
they can be treated as a single anyon. The density matrix of the target system is 



a,a' ,c,c' ,f,fi,^' 



P{a,c-J,t,){a',c'-jy) „TM 
/ „xl/4 ■'Tm' 



a,a' ,c,c' ,f,fi,fi' 



dada' dcdc' (PfJ a' 




(3.3) 



We will assume that the probe anyons are also not entangled with each other, and 
that they are all identical (or, more accurately, belong to an ensemble of particles 
all described by the density matrix p^). We will consider generalizations of the 
probe anyons in Section [331 Such generalizations complicate the bookkeeping of 
the calculation, but will have qualitatively similar results. A probe system involves 
the probe anyon B, which is sent through the interferometer entering the horizontal 
leg s and the anyon D which is entangled with B and will be sent off to the 
(left) side. We will write the directional index s of the probe particle as a subscript 
on its anyonic charge label, i.e. bg. The density matrix of a probe system is 
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= E pfd,b^;h,x){d',bU;h,x')^ h, A) {d', /i, A'| 



b,b' 44' ,h,\,X' 



y-v P{d,b^;h,X){d',b'_;h,X') ^Ya 

b,b' 44' ,h,x,x' (dddd'dbdb'dl)^^^ d'x\b' 



(3.4) 



The unitary operator representing a probe anyon passing through the interferometer 
is given by 

U = T2ET1 



e^'-R-/^ 



c^'^Rba 

This can be written diagrammatically as 



(3.5) 
(3.6) 



U 



J6i 



—f-,f* —r*f* 
LIL2 I il^2 



B 



\ 



+e 



i0u 



A 



rit2 -tlt2 



V 



(3.7) 



The position of the anyon C with respect to the other anyons must be specified, and 
we will take it to be located below the central interferometry region and slightly to 
the right of A. (The specification "slightly to the right" merely indicates how the 
diagrams are to be drawn, and has no physical consequence.) For this choice of 
positioning, the operator 



V 



RcB 
RcB 



y 

B ^\ 



(3.8) 



represents the braiding of C with the probe. In Section [J4l we will give the results 
for situating the anyon(s) C in different locations outside the central interferometry 
region, and find that the qualitative behavior is essentially the same. 

After a probe anyon B is measured at one of the detectors, it no longer interests 
us, and we remove it along with its entangled partner D from the vicinity of the 
target anyon system. Mathematically, this means we take the tensor product of the 
probe and target systems, evolve them with VU (which sends the probe through the 
interferometer) to get 

p = rf/ (p^ ® p^) f/v^ (3.9) 

apply the usual orthogonal measurement collapse projection 



Pr(s) =Tr[pn, 
1 



P 



Ft (s 



-n,pn. 



(3.10) 
(3.11) 
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with = \s) {s\ for the outcome s, and then finally trace out the anyons B and D. 
Since the probe anyons are all initially unentangled, we may obtain their effect on 
the target system by considering that of each probe individually. 



3. 1 One Probe 



We begin by considering the effect of a single probe with definite anyonic charge h. 



I.e. p" 



6, h^-l 



, and return to general immediately afterwards. For 
a particular component of the target anyons' density matrix, the relevant diagram 
that must be evaluated for a single probe measurement is 



6" 



a 



u 



(3.12) 



For the outcome s this is 
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1^ 



u 



(/,/.,/.')(e,",/3) 



U 
an 



a' d 

e,a,l3 



w6_ 



X < 



Itil'hl' 



a c 

n 



a c 



CD 



a' d 



a c a c 

I g I 



4E [(^av)"' 



a' c' 



(f,li,li'){e,a,f3)' 



A e A 



a c 



♦ e ♦ 



a c 



e,a,/9 



a' c' 



where we have defined 



P:a'e,b = \tif \r2f + tir^r^t^e^^^'-^-^M^b 

+ttrit2r2e-'(^'-^")M:,, + |ri|' l^sl' (3.14) 

and have used Eqs. (12.3712.381) to remove the b loops. A similar calculation for the 
s =T outcome gives 



-tlnt2r2e-'^''-'''^M:,, + \nf Ir^f . (3.15) 

From this, inserting the appropriate coefficients and normalization factors, we find 
the reduced density matrix of the target anyons after a single probe measurement 
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with outcome s: 



Pr (s) 



EP{a,c;f,fM){a',c'-,f,i-i') Paa'e,b 
( 2\ Pr (s) 

a,a',c,c',f,fj,,ii' [dada'dcdc'd f) ^ ' 



X 



E 



aa'e.b 




, {dfdff^ Pr(s) 
X K'c'](e,a,/3)(/',.,.') c; i/) (a', c'; i/'l 



a,a' ,c,c' ,f ,11,11' 
e,a,(3,f',v,u' 



(3.16) 



where the probabiUty of measurement outcome s is found by additionally taking 
the quantum trace of the target system, which projects onto the e = 1 components, 
giving _ 

Pr (s) = Tr [pn,] = pta,c-,M(a,c-,MPLi,b- (3- H) 



a,c,f,fi 



We note that 



Prai,6 = l^il' \r2\' + \rif \t2f + 2Re {hrlr;t;e'^'^-'-^ M^b} (3.18) 
l^il' l^2|' + \rif \r2f - 2Re{tirtr;t;e*('^'-^°)M„b} (3.19) 

give a well-defined probability distribution in s (i.e. < Pa^i < 1 and Paai,b + 
The quantity 

tirlt*r*e'^^'-^''^ = Te'^ (3.20) 

determines the visibility of quantum interference in this experiment, where vary- 
ing 9 allows one to observe the interference term modulation. The amplitude T — 

\tirit2r2\ is maximized by \tj\ = \rj\ = l/\/2. In realistic experiments, the ex- 
perimental parameters tj, rj, 9i, and On will have some variance, even for a single 
probe, that gives rise to some degree of phase incoherence. Averaging over some 
distribution in 9, one finds that e*^ in the interference terms should effectively be 
replaced by (|e*^^ = Qe*^*. In this expression, e'^* is the resulting effective phase, 
and Q e [0, 1] is a suppression factor that reflects the interferometer's lack of coher- 
ence, and reduces the visibility of quantum interference. For the rest of the paper, 
we will ignore this issue and assume Q = 1, but it should always be kept in mind 
that success of any interferometry experiment is crucially dependent on Q being 
made as large as possible. 
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We can now obtain the result for general by simply replacing p^^'e b everywhere 
with 



faa'e,B = B (&) PL'e,b 

b 

P^sib) = ^ P{d,b^;h,\){d,b->;h,X)- 
d,h,X 



(3.21) 
(3.22) 



We will also use the notation Mas = Jlb^^B (b) Mab- That this replacement gives 
the appropriate results follows from the fact that we trace out the D anyon, and may 
be seen from 



Tr 



D 



P{d,b^;h,X)(d,b'_^;h,X') 

b,b' ,d,h,X,X' {d'^dbdh'dl)^^^ 




Pfd,b-,;h. 



A)(d,f)_;/i,A)^ 



h 6_ 



P{d,b^;h,^){d,b^;h,X) ^ 



bA,h,X 



b 6_ 



■E^rn (&)Tr5 

b 

:Tr^^Pr5(6) 



b,b^;l)(b,b^;l 



b,b^;l)(b,b^;l 



(3.23) 



where we used Eq. (12.271) in the first step. 



3.2 N Probes 



The result for initially unentangled identical probe particles sent through the in- 
terferometer may now be easily produced by iterating the single probe calculation. 
The string of measurement outcomes (si, . . . , s^) occurs with probability 



Pr (Si, . . . , Sn) — ^ P{a,c;f,ti){a,c-J,ii)Plal,B ■ ■ ■ Paal,B 
a,cj,fj. 



(3.24) 



and results in the measured target anyon reduced density matrix 
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P (Si,..., Sat J 



P{a,c-J,^i){a',c'-J,^i') Paa'e,B ' • ■ Paa'e,B 



E 



(/,M,M')(e."./3) 



K:'](e,«,/3)(/',.,.') h c; /', u) (a', c'; /', u'\ . (3.25) 



It is apparent that the specific order of the measurement outcomes is not important 
in the result, but that only the total number of outcomes of each type matters, hence 
leading to a binomial distribution. We denote the total number of Sj =^ in the 
string of measurement outcomes as n, and cluster together all results with the same 
n. Defining (for arbitrary p and q) 



Wn (n; p, q) 



N\ 



pnqN-n 



n\ {N - n)\ 

the probability of measuring n of the N probes at the horizontal detector is 



Pr^(n)= Y: 

a,cj,fj. 



(3.26) 



(3.27) 



and these measurements produce the target anyon reduced density matrix 



pi (n) 



E 



Pi 



I t £ I 

a,a ,c,c ,j 
e,a,f3J',uy 



{dfdf,f'^ Pr^ in) 

c- /', v) (a', c'; /', v'\ . (3.28) 



In Ref. [[82ll . we obtained the reduced density matrix that ignores the measurement 
outcomes and describes the decoherence (rather than the precise details of collapse) 
due to the probe measurements. We find this density matrix by averaging over n, 
giving us the result in Eq. (15c) of [|82J, though for more general target and probe 
systems 



N 



pi = T. Pi^^ (^) pi (^) 

n=0 



E 

a,a',c,c',f,ii,fi' 

-1 



tif MeB + \rif) 



N 



where we used 



P{a,c;Mia',c';fy) 
[dfdrfl' 



(3.29) 
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W^AT {n',Paa'e,B^Pla'e,B) ~ (Paa'e,B + PL'e,^) 

n=0 

= (|ti|'M,B + |ri|')'^. (3.30) 

The interferometry experiment distinguishes anyonic charges in the target by their 
values of p^^i ^, which determine the possible measurement distributions. Differ- 
ent anyonic charges with the same probability distributions of probe outcomes are 
indistinguishable by such probes, and so should be grouped together into distin- 
guishable subsets. We define for /t = 1, . . . , m < \C\ to be the maximal disjoint 
subsets of C such that p^^^^ b = for a &C^, i.e. 

C^ = [aeC: p-2ai,B = P^} ■ (3-31) 
Note that p^ai b = Pa'a'i B (for two different charges a and a') iff 

Re [tirlrltle'^^'-^''^ MaB) = Re {hrlr^t^e'^^'-^''^ Ma'B} (3.32) 
which occurs either when: 

(i) at least one of ti,t2,ri, or r2 is zero, or 

(ii) \MaB\ COS {d + Lpa) = \Ma'B\ COS {9 + ifa'), whcrc e = arg (tir^fr^t^e^^^'"^")) 
and ifa = arg (M^s). 

If condition (i) is satisfied, then there is no interference and Ci = C (all target 
anyonic c harg es give the same probe measurement distribution). Condition (ii) is 
genericalljiil only satisfied when MaB = Ma'B, but may be non-generically satis- 
fied by setting 9 = — arg {MaB — Ma'B} ± f • With this notation, we may write 

Pr;v (n) = ^ PiA (n) Wn (n; 1 - p^) (3.33) 

K 

Pta{k)= Pta,c;M{a,c;M- (3-34) 

a£CK,c,f,fJ. 

We emphasize that if the parameters tj, rj and 9 in the experiment are known and 
adjustable, then the measurements may be used to gather information regarding the 
quantities Mat, which, through its relation to the topological S'-matrix, may be used 
to properly identify the anyon model that describes an unknown system [185]| . 

In Sections [3.2. ll and B. 2. 2[ we will show that, as ^ oo, the fraction r = n/N of 
measurement outcomes will be found to go to r = with probability Pr^ (k), and 

^ The term "generic" is used in this paper only in reference to the collection of interferom- 
eter parameters tj, rj, 0i, and On- 
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the target anyon density matrix will generically collapse onto the corresponding 
"fixed states" given by 



Pk = 2^ , ,1/2 ^aa'e.B [Pk] 

a,a',c,c',f,tM,tM' {".fd-f) 



X 



e,a,l3J',uy 

U,.y)ie,a,,) Kc'](e,«./3)(/',.,.0 ^5 -) ^'1 (3-35) 



where 



. / A ; P?fe) = ^ ~ P--'-'B = and a, a' G 

Aaa'e,B (Pk) = <( ^ . (3.36) 

otherwise 

Fixed state density matrices are left unchanged by probe measurements. We also 
emphasize that the condition: p^a'e s = ^ ~ PL'e b = is equivalent to M^b = 
1 (noting that Mg^ = 1 implies MaB = Ma'B and a, a' G C^). This gives the 
interpretation that the probes have the effect of collapsing superpositions of anyonic 
charges a and a' in the target that they can distinguish by monodromy (MaB 7^ 
Ma' b), and decohering all anyonic charge entanglement between the target anyon A 
and the anyons C outside the central interferometry region that the probes can "see" 
by monodromy, i.e. removing the components of the density matrix corresponding 
to e-channels with M^b 7^ 1- Non-generically, it is also possible to collapse onto 
"rogue states," for which the diagonal density matrix elements are all fixed and 
some of the off-diagonal elements have fixed magnitude, but phases that change 
depending on the measurement outcome (i.e. are "quasi-fixed"). Because rogue 
states occur only for specific, exactly precise experimental parameters, they will 
not actually survive measurement in realistic experiments. We note that if M^b = 1 
only for e = 1, then the probe distinguishes all charges, and the fixed states are 
given by 

A >r-Pr^(c|a) ^ ^ PrA(c|a) , , 



where 



Fta (c|a) = ^^-^ , (3.38) 

P{a,c;f,f,)ia,c;M 

c,f,tJ. 

for which the target anyon A has definite charge and no entanglement with C. We 
give examples of fixed state density matrices for several significant anyon models 
in Section [3]6l 

In principle, one may also consider the "many-to-many" experiment described in 
Ref- L83J5 where the target anyon system is replaced with a fresh one (described 
by the same initial density matrix) after each probe measurement. For this type 



26 



of experiment, the result for each probe is described by the single probe outcome 
probability, Eq. (13.171) : 

^^^(■5)= X! Pta,c;f,fi)(a,c;f,fi)Plal,B- (3.39) 
a,cj,fj. 

Thus, for such probe measurements, the number n of s =^ measurement out- 
comes will have the binomial distribution: Wn (n; Pr (^) , Pr (t)). 



3.2.1 Large N 

We would like to analyze the large N behavior of the measurements. This is essen- 

tially determined by Wn {n; p^, 1 — p^) and — Pri^'(n) ' ^'^ ^'^^ 

sider these in detail. Of course, Wn (n; p^, 1 — p^) is just a familiar binomial dis- 
tribution. Changing variables to the fraction r = n/N of total probe measurement 
outcomes in the horizontal detector, the distribution in r is given by 

WN (r; p„ l-Pn) = Wn {rN- p,, l-p^)N (3.40) 

and has mean and standard deviation 

r=p. (3.41) 

Ar = a^ = ^p^{l~p^)/N. (3.42) 

Taking large and using Stirling's formula, this may be approximated by a Gaus- 
sian distribution 

wtv (r;p«,, 1 -p«) =i ^^=e ^"'i- . (3.43) 



Taking the limit ^ oo gives 

lim Wn (r; p^, I - p^) = S {r - p^) (3.44) 

(defined such that 5 [r — p) dr = 1, when p = or 1), so the resulting probability 
distribution for the measurement outcomes is 

Pr (r) = lim Pr^v (r) = V Pr^ (/?) 5 {r - p^) (3.45) 

Thus, as ^ oo, we will find the fraction of measurement outcomes r ^ p^^ with 
probability Pr^ (k). 

Though the probability of obtaining the outcome r which is away from the closest 
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vanishes as 

^^A. (r; l-p^)^^/Ni ^^i,^ (3.46) 

\ [1 — r) ) 

for large A^, the resulting density matrix should still be well defined for all r (at least 
for large, but finite A^). In particular, we will use the positivity property of density 
matrices, in the form of the Cauchy-Schwarz type inequality p^^iPw > \p^ii^\ , to 
evince their large N behavior in terms of conditions on p^^'e b- From the quantity 



^N:aa'e,B (r) = 



Wn [rN;p^^, B,Paa'e,B 



FiN (rN) 



Pk' 



\Paa'e,B , 



1 - Pk' 
, Paa'e,B 



1-r 



(3.47) 



we can see that as ^ oo, the e = 1 terms (those that determine the "diagonal" 
elements) behave as: 



(i) AN;aai,B (r) p^jij^ « G C^, , if Pi A 7^ and p^^ (1 > 
{l-p^y^'' for all K ^ Ki, 



(ii) Aiv-aai R (r) 



1 



for aeC^.U C^,, if Via (/ti) + Pr^ {^^2) 7^ 



andp;Jl-p«J' '^=p;Jl-p^J^ '^>pr(^_pji '■forallK ^ Ki,K2, or 
(iii) AAT-aaiB ^ for a G C^i, if there is some k with Pr^ (k) ^ and 

(1 - vk) > pi;, (1 - P^^f^''- 

If a G C^i, where p;:^ (1 - Pn^f^'' > pl{l- Pk)^''' for all k ^ ki, but Pr^ (ki) = 
0, then A7V;aai,B ('^) ^ cxd. However, Pr^ (ki) = also implies that the density 
matrix coefficients involving a are strictly zero, so we need not worry about this 
case. 



We note that for each k, the variable r has a closed interval J^, containing in its 
interior, such thatp^ (1 — Pk)^ > p^/ (1 — Pk')^ ' n' k (i.e. /«; satisfies 

(i) in its interior and (ii) at its endpoints). We say that r is congruous with in this 
interval /«; (and congruous to two different at the intersecting endpoints of such 
intervals). 

For arbitrary (in particular, the "off-diagonal") terms, the positivity condition com- 
bined with Eq. (13.471) as ^ oo tells us that we must have /^N;aa'e,B {r) 0, 
except when r is congruous with both a and a' , in which case \/S.iq.^aa'e,B < 
^N;aai,B (r), and Aiy-aa'e,B (^) oo should uot bc allowcd (except when the den- 
sity matrix elements involving a or a' are strictly zero, making it irrelevant). From 
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this we find that either: 



(a) there is some n (possibly even with a and/or a' in C^) with Pta (k) 7^ and 



Paa'e,B Pla'e,B ^ ' in which CaSe AN;aa'e,B (r) ^ 0, Or 
~ {Paal,B^ (pLi.s) — (Pa'a'l,s) {Pa'a'l,B^ 

r congruous with both a and a' , in which case | AjV;aa'e,B {t) \ — ^N;aai,B {t) 
Case (b) deserves some further inspection. First, we note that we have 



1-r 



(b) 



Paa'e,B 



P 



T 

aa'e.B 



1-r- 



Paa'e,B 



T 

Paa'e,B 



l-r 



<p:(i-p.y 



(3.48) 



on the entire interval congruous with a e C^. If there is some point r* in the 
interior of h for which 



— ^ 


r* 


T 


1— r* 


Paa'e,B 




Paa'e,B 





1— r» 



(3.49) 



then in order not to violate the inequality when r is increased or decreased from r*, 
we must have 

Paa'e,B 



It follows that 



Paa'e,B 



Paa'e,B Paa'e,B 



Pk 
1 



1-r 



1-r 



on the entire interval /„, and, more significantly, that 



Paa'e,B 



1 - 



Paa'e,B 



Pk- 



(3.50) 
(3.51) 
(3.52) 



The same argument holds with respect to a' instead of a, giving the additional 
condition a, a' G C^. Hence, even at exponentially suppressed r, superpositions of 
anyonic charges from different do not survive measurement. 



Pushing this a bit further, we note that for r e [0, 1] and fixed p e [0, 1] 



r"- {I - r^-'' > p"- {1 - p) 



1-r 



(3.53) 



with equality atr — p. The positivity condition gave us (rewriting (a) and (b)) 



Paa'e,B 



with equality for r e if a, a' e and Paa'e,B 
these, we have 

1 ^ 

r'' (1 - r) > 
for all r, with equality occurring at r = p^„,g ^ 



T 

Paa'e,B 
1- 



1-r 



T 

Paa'e,B 



(3.54) 
Pk. Combining 



T 

Paa'e,B 



1-r 



T 

Paa'e,B 



. If 



Paa'e,B 



+ 



1 

Paa'e,B 



only when Pa ^p^/ 
> 1, then there is some r (e.g. r = 



(3.55) 

Paa'e,B 
Paa'e,B 



) 
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for which r^' (1 — r 



l-r 



< 



T 

Paa'e,B 

we have 



1, then r*" (1 — r 



Paa'e,B 
l-r 



l-r 



Paa'e,B 
Paa'e,B 



, violating Eq. (13.551) . If 

l-r 



Paa'e,B 



+ 



Paa'e,B 



at r 



Paa'e,B 



Paa'e,B 



+ 



Paa'e,B 



< 1 



. Hence, 



(3.56) 



with equality only if pa = Pa' 



Paa'e,B 



Paa'e,B 



. We believe one should 



be able to to show that this condition on p^^^'e b follows directly from the properties 
of anyon models, in which case these arguments could be made in the opposite 
direction, i.e. that positivity of the density matrix being preserved by these probe 
measurements follows from properties of anyon models; however, we have been 
unable to succeed in doing so. 

For P2a'e,B = PacC'^^-'-^.s and pL'e,B = (1 - Pfc) e'^--'--B, we see that if aaa'e,B = 

(3aa'e,B, thcU 



\h\'M,B + \ri\=P:a'e,B+Pla'e,B = e*"""'='" 



(3.57) 



implies that either: (a) ri = and M^b = e*"'''''^ ^ , or (b) M^s = 1 and aaa'e,B = 

f3aa'e,B = 0. 

One might also find it instructive to consider a large expansion (using Stirling's 
formula) around p,^ to get 



(r; Paa'e,B, Pla'e,B) - {r; p^, 1 - p^) 6 



-Gm ( r;p^ , , 

V aa'e,B aa' e 



^N:aa'e,B (r) 



1 -Gn 

-e 



i T"X) * X)^ 



(3.58) 
(3.59) 



Pk 



Gn {r;p,q)^N 

+N (r — p^) In 



Pf,\n l—j + (1 -p^) In 



1 -Pk 



, P (1 -Pk) 



(3.60) 



Clearly, e '^Nir;p,g) gives exponential suppression in A^, unless p = PkC^" and q 
(1 — Pk) e*^, in which case 



g-Gjv(r;p,g) ^ ^i[ap^+l3{l-p^)]N ^i{a-l3)N{r-p^) 



(3.61) 



(which is equal to 1, when a = P = 0). We also note that integrating the quan- 
tity Wat (r; Paa'e,B^ Pla'e,B) o^^^" vanishcs exponentially in N, unless aaa'e,B = 
Paa'e,B = 0, wMch is why such quasi-fixed terms do not appear in Eq. (13.291 ). the 
density matrix obtained by ignoring measurement outcomes (except in the case 
when ri = 0). 

To summarize, we found that, for large A^, the quantity /^N;aa'e,B {r) vanishes ex- 
ponentially unless r is congruous with a, a' G and p^a'e b 



T 

Paa'e.B 



Pk- 
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This means a measurement outcome fraction r exponentially collapses the density 
matrix onto one that has support only in C^, and consequently will drive r toward 
Pk. The resulting target anyon reduced density matrix 



E 



P(a,c;/,M)(a',c';/,/.') 



aa'e,B 



(r) 



a,a' ,c,c' ,f,fj.,fj,' 



(/,M,M')(e,",/3) 



where 



Aaa'e R (r) 



lim AAr.<j„/e,B (r) 



is found with the probability distribution 

Pr (r) = ^ Via 5 (r - . 

K 

The resulting density matrices are of two forms: 



(3.62) 
(3.63) 

(3.64) 



{I) fixed states, for which all non-zero elements of the density matrix correspond 

to V2a'e,B = 1 - Vla>e,B = Vn, and 

(2) rogues states (or quasi-fixed states), for which all elements of the density matrix 



correspond to 



Paa'e,B 



Paa'e,B 



Pfi, but for some of the "off-diagonal" 
elements (e ^ 1) withp;^^,^^^ = pi^e^^'^"'^-'^ ^^'^Pla'e,B = (1 ~" Pfc) e*^""'"'^' where 

Oiaa'e,B ^ud (3aa'e,B ^rC UOU-ZCrO (UUlcSS Ti = 0). 

Fixed states have the property that probe measurements leave their density ma- 
trix invariant. Rogue states have the property that probe measurements leave their 
"diagonal" elements and possibly some of their "off-diagonal" elements invariant, 
while some of their "off-diagonal" elements are unchanged in magnitude, but have 
a changing phase. We will see in Section 13.2.21 that satisfying the conditions for 
rogue states requires non-generic experimental parameters. 



3. 2. 2 Minding our p's 

Li Section 13.2.11 we have shown that performing many probe measurements col- 
lapses the target density matrix onto its elements which correspond to pia'e,B satis- 
fying 

Paa'e,B\ = 1 " \pL'e,B\ = Pk (3.65) 

for a, a' G C^, so we would like to determine when this condition is satisfied. 



For completeness, we first list the results for the trivial cases where there is no 
actual interferometry (for which Ci = C): 
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(i) When ti = 0, we have p2a'e,B = I'^sl^ and p\a'e,B = ksl^, so all elements are 
fixed. 

(ii) When ri = 0, we havep^^/g ^ = \r2\^ M^b andp^^,g ^ = \t2f M^b, so elements 
with M^B = 6*'^'=^ {(feB 7^ 0) are quasi-fixed, and those with M^b = 1 are fixed. 

(iii) When ^2 = (and h ^ 0), we have Paa'e,B = ^eB and pla'e,B = 

so elements with MgB = e**^"^^ i'feB 7^ 0) are quasi-fixed, and those with MeB = 1 
are fixed. 

(iv) When r2 = (and ti ^ 0), we have p'^a>e,B = l^if and pL'e.s = l^il^ ^eB, 
so elements with MeB = e*'^'^^ {'feB 7^ 0) are quasi-fixed, and those with M^b = 1 
are fixed. 

From here on, we assume that |tirit2?"2| 7^ (unless explicitly stated otherwise). 
We begin by considering the more stringent condition necessary for fixed elements. 
Using p2a'e,b + Pla'e,b = hl^ ^eB + and Eq. (12.401) . we have: 

(v) (When ti ^ 0) An element is fixed, withp^^,^^ = l-p\a'e,b = P^, iff MeB = 1, 
and this implies MaB = Ma'B and a, a' E C^. 

Thus, even without initially requiring a, a' E (from positivity), we find that it is 
a necessary condition for fixed elements. 

Now, we examine the conditions that give quasi-fixed elements. Such terms have 

Paa'e,B = P^e'" and p^^/^^g = (1 e*'^''-'<=.s, with aaa'e,B 7^ (3aa'e,B and 

a, a' E C^. (Recall that if aaa'e,B = f3aa'e,B and ri ^ 0, then aaa'e,B = f3aa'e,B = 0.) 

Examining these conditions for Mas = Ma'B, we find 



0=|t2 



Paa'e.B 



2 \ I I |2 

Pk + m 



p 



T 

aa'e,B 



Pk) 



1+ |4 \± |2 I |2 
\h\ ^2 



(|MeB| 



+2 ItiHnl' |t2|' |r2|' (Re {MeB} - i; 



(3.66) 



which requires MeB = 1 and, hence, gives us: 

(vi) (When \tirit2r2\ 7^ 0) There are no quasi-fixed elements for Paa'eB with 
MaB = Ma'B, Only fixed ones. (In particular, this applies to a = a' .) 

For MaB 7^ Ma'B, wc Can only have a, a' E (i.e. plai.B = Pl'a'i,B) when the 
experimental parameters are tuned to = — arg {MaB — Ma/^} ± f , so quasi-fixed 
elements only occur non-generically. From the conditions on Paa'e B' these values 
of ^, we find 
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= P'2a'e,B -pI- Pla'e,B + (l " 

-2 (1 - Re {MeB}) 2 \hnt2r2\Rc {e^'MaB} 
+2 \ti\^hn{MeB} \hnt2r2\hn{e'''MaB + e-'^M^r^} 



(3.67) 



and 



0= \h 



Paa'e,B 



2 \ , I |2 



Paa'e,B 



'^-Pk] 



= \tit (iMesl' - 1) + 2 Inl' 1^2!' |r2|' (Re{MeB} - 1) 

+ {\tint2r2\ Im {e'^MaB + e-^^M„%})' (3.68) 

which may be rewritten to give: 

(vii) Quasi-fixed elements with ^ only occur non-generically, and the condi- 
tions (when \tirit2r2\ 7^ 0) that must be satisfied for them to occur are: 



^ = -arg{M„5-M„,s}± 2 



(3.69) 



Im {e''MaB + e-^'M^s}]' = (l - \M,b\^) + 2 (1 - Re {M,b}) (3.70) 



Re{e*^M„B} = 



4 IrJ 



r4-f4 (1 



r2\ 
t2\ 



eB\ 



+ -Im {MeB} hn {e'^MaB + e-'^M^^}] (1 - Re {Mes})"' . (3.71) 



To demonstrate that it is, in fact, sometimes possible to satisfy the conditions for 
quasi-fixed elements given in (vii), we present the following example: 

Consider an anyon model which has at least two different Abelian anyons a and 
a', and some anyon b for which Mab = e*'^"'' and Ma'b = e^'^'^'^ are not equal (for 
example, almost any Z^r model, such as Z2^^^^ or Zg^^ is sufficient). The difference 
charge e is uniquely determined (since a and a' are Abelian) and has Mg^ = e*'^^'' = 
eii.fab-fa'b) , Setting 9 ^ -\ {cpab + (fa'b) + mi gives 
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aa'e.b ' 



T 

Paa'e.b ' 



ti\ \r2\ e 



+ \ri\ \t2 



-\ti\\t2\e'y^^-") + \n\\r2 



+n7r) 



Paa\,b —Pa'a'l,b 



Paa' e,b 



Pa 



\tif' ^2!^ + 2 \tirit2r2\ cos + /itt) + |ri|^ 1^2!^ 



(3.72) 
(3.73) 

(3.74) 



In fact, it turns out this example is the only way to satisfy the conditions for quasi- 
fixed elements with [Mesl — 1. Indeed, this can even be shown without initially 
requiring a,a' E from positivity. It seems rather difficult to satisfy the conditions 
for quasi-fixed elements when \MeB\ 7^ 1, and we suspect (but are unable to prove) 
that it may, in general, actually be impossible. It is certainly not possible to have 
quasi-fixed elements with jMesI 7^ 1 for arbitrary non-Abelian anyon models, as 
one can check that they do not exist for either the Ising or Fib anyon models, for 
example. 



3.3 Distinguishability 



We would like to know how many probe anyons should be used to establish a de- 
sired level of confidence in distinguishing between the various possible outcomes. 
For a confidence level 1 — a, the margin of error around is specified as 



(3.75) 



i.e. the interval [p^ — Efi,pi^ + contains 1 — a of the probability distribution, 
where 2;* ^2 defined by 



1 — a 



(3.76) 



To achieve this level of confidence in distinguishing two values, pi and p2, we pick 
N so that these intervals have no overlap 



Ap = |pi - P2I >Ei + E2 = z*^i2 (0-1 + 0-2 

— y* 

— ^al2 



N 



N 



which gives the estimated number of probes needed as 



N > 



^a/2 + ^JP2{^-P2)) 



Ap 



(3.78) 
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Since p{l — p) < |, we could conservatively estimate this for arbitrary pj as 



N > 



(3.79) 



1^ ~ 1^2!^ ~ <S 1, and Ap is of 



On the other hand, if pi and p2 are of order \t^ 
order 2t'^AM, where AM = \Ma^B — Ma^s], (i-e. employing 9 such that Ap is as 
large as it can be) then we can estimate 



~ \tAM J 



(3.80) 



We note that for any two outcome probabilities, pi and p2, there are always two 
values of 9 (i.e. non-generic conditions) that make pi = p2, and hence indistin- 
guishable. Here are the values of z^^^^ some typical levels of confidence 



1 — a 


.6827 


.9545 


.99 


.999 


.9999 




1 


2 


2.576 


3.2905 


3.89059 



For greater confidence, the number of probes needed roughly scales asN — log a. 

A special case of interest exists when = 1 and jAfaj^l = 1 for one of two 
probabilities that we wish to distinguish. In this case, using 9 = tt — arg{Ma^s} 
gives pi = 0, so any measurement outcome s =— > automatically tells us the target's 
anyonic charge is not in Ci. If the alternative outcome has p2 7^ 0, 1, then Ci and 
C2 are said to be sometimes perfectly distinguishable, since a s outcomes 
tells us the target's anyonic charge is in C2. If M^^b = —Ma^B and we also have 
\tjf = 1/2, then p2 = I, and Ci and C2 are always perfectly distinguishable, since 
any single probe measurement will indicate whether the target's anyonic charge is 
in Ci or in €2- 



3.4 Target System Configuration 



In this section, we consider the effect of locating the anyons C that are entangled 
with the target A in different regions outside the central interferometry region. If C 
is located above the central interferometer region, we would have 



V 



Rbc 
Rbc 



(3.81) 



for which similar diagrammatic evaluation gives 
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P:a'e,B = \ti\' \r2\' + hrlvltle'^'^-'-^ M,, B 



(3.82) 



(3.83) 



instead of Eqs. (|3. 1413. 1513.211) . If C is located between the output legs of the 
interferometer, we would instead have 



V 



Rbc 





(3.84) 



The resulting diagrammatic evaluation in this case gives 



+ttrit2r2e-^(^'-^")M:B + Itsi' M:^, (3.85) 
p\a'e,B = \tif \t2f M,b ' hrlrltle^^'^-'^^ M^b 

-t*rit2r2e-'(''-'")M:,^ + |ri|' jral' . (3.86) 

For both of these cases, the arguments from before apply directly and limiting be- 
havior is exactly the same. One can also envision more complicated situations, such 
as having the C anyons distributed amongst all the regions outside the central one. 
The resulting calculations are straightforward, but too cumbersome to display here 
explicitly. However, the limiting behavior is essentially the same as before, as one 
would expect from the previous analysis: Interferometry measurement generically 
collapses the target system onto fixed states, which are characterized as having the 
target anyons A (those which the probes interfere around) in a charge subset that 
the probe cannot distinguish by monodromy (i.e. a, a' such that MaB = Ma'b), 
and the anyons in distinct regions having no coherent anyonic charge entanglement 
crossing the probe anyons' beam paths that the probes can "see" by monodromy 
(i.e. anyonic charge entanglement characterized by charge e can entangle anyons in 
distinct regions only if MeB = !)• 



3.5 Probe Generalizations 



In this section, we examine the effects of using probe systems that are even more 
general than those employed so far. We will first consider generalizing the input 
direction, so that probes may enter in arbitrary superpositions of the two input di- 
rections. Then we will consider the use of probes that are not identical, so that each 
probe system is described by a different density matrix. For both of these, the probe 
systems and target system are all still initially unentangled. One may also consider 
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cases where there is nontrivial initial entanglement between these systems, or post- 
interferometer charge projections, but these typically lead to qualitatively different 
behavior, and greatly increase the complexity of analysis, so we will not consider 
them here. 

3. 5. 1 Generalized Input Directions 

For probes that are allowed to enter the interferometer through either of the input 
legs, possibly even in superposition, the probe systems' density matrices take the 
form 

b,b' ,d,d' ,h,\,X' ,r,r' ^ ' "-h 

Using this, we find the same result as before, except, instead of Eqs. (13. 1413. 1513.211) . 
the values of pla'e,B given by 

Paa'e,B — X! P{d,br;h,X){d,b^, ■,h,^)Paa' e,b,r,r' (3.88) 
d,h,X,b,r,r' 

where 

P:a'e,b,^,^ = Ihf \r2\' + t^rltlr^e^^'^-'-^ M^t 

+tlnt,r,e''^'^''''^M:,, + \nf \t,f (3.89) 
PweA-,! = ^1^1 l^2|' Meb - hhtlr^e'^'^-'^'^ M^b 

+nnt2r2e-'^''-'''^M*, - tiri Ital' (3.90) 
Pl'eM,^ = tlrl Ir^f M,, + rtrtt;r;e^(^'-^")M,, 

-tltlt2r2e-'^'^-'''^M:,, - tlrl 1' (3.91) 
P:a'eAt,t = Me, - t^rlt^r^e'^'^-'-^ M^t 

-tlnt2r2e-'^''-'''^M:,, + \t2\' (3.92) 

and 

Pla'e,b,^,^ = \tif M,, - tirtt;r;e^(^'-^")M,, 

-ttrit2r2e"'(''-'")M:,, + \nf |r2|' (3.93) 

-ririt2r2e~'(^'~^")M*f, - tin |r2|' (3.94) 
Pla'eM,^ = tVi \t2f Me, - rlrHlr^e'^'^-'-^ M^b 

+fifit2r2e-*(''-'"^M:,, - tlrl |r2|' (3.95) 
Pl'eAU = l^il' Me, + hrltlr^e'^'^-'-^Mab 

+tlnt2r2e-'^'^~'''^M:,, + \t,\' \r2f . (3.96) 
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It is straightforward to check that 

Paal,B Plal,B = X! Pfd,br;h,X){d,br;h,>d ~ (3.97) 

d,h,\,b,r 

and one can see that, generically, the only terms in the target anyons' density matrix 
that will survive many probe measurements are those in e-channels with 

MeB= E pfdA;h,mA;h,x}^eb = l. (3.98) 



3.5.2 Non-Identical Probes 

When the probes Bi, . . . , are described by different density matrices p^^ (though 
are all still unentangled with each other and with the target system), we must use 



PL'e,B,-j:P'^BAb)pla'e,b (3-99) 
b 

Prs, W = E p5V;M){<i.6.;M) (3. 100) 

d,h,X 

for each probe. This gives us the probability for the string of measurement out- 
comes (si, . . . , Sat) to occur as 

Pr (Si, . . . , Sn) = E ^(l,c;/,|u)(a,c;/,/i)Paal,Bi • • -Paal.B^v' (3.101) 
a,c,f,iJ. 

with the resulting target anyon density matrix 



a,a' ,c,c' ,f,fMy 



y{a,c■J,^l){a',c'■,f,^l') Paa'e,Bi ■ ' ■ Paa'e,BN 



(/,/^,/i')(e,a,/3) 



{dfdf^f^ Pr(si,...,sjv) 
Kv](e,a,/3)(/',.,.') /, v) {a', c'; . (3.102) 



With this generalization, we find that the order of measurement outcomes does, in 
fact, matter. This is obstructive to providing a quantitative description of the large N 

behavior; however, the qualitative behavior should be transparent after the analysis 
in previous sections for the identical probes. Each probe measurement will execute 
some amount of projection, to some extent collapsing superpositions of anyonic 
charges that the probe is able to distinguish by monodromy. 
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3.6 Examples 



In this section, we apply the general results to some important examples, specif- 
ically: the Zjv, Fibonacci, and Ising anyon models. (The application to some ad- 
ditional important examples, such as SU (2)^ and D (Z^v), may be found in 1791 .) 
All of these have N^i^ = 0, 1, so we will drop the fusion/splitting spaces' basis la- 
bels (greek indices), with the understanding that any symbol involving a prohibited 
fusion vertex is set to zero. Anyon models are completely specified by their F- 
symbols and i?-symbols, so we will provide these, as well as list some additional 
important quantities that can be derived from them, for convenience. To relate these 
to interferometry experiments, we give the corresponding fixed state probabilities 
Pk and density matrices p^, as described in Section [321 

3.6.1 zj;) 

The Abelian Zat anyon models ^ have anyonic charges C = {0, 1, . . . , — 1}, 
where here designates the vacuum charge. The fusion rules are just given by Zat 
addition, and, to denote this, we define [n]j^ E C as the least residue of n mod A^. 
The Z^'' anyon models, where w = n for A^ odd and w = n and n + | for A^ even, 
with n = 0, 1, . . . , A^ — 1, are described by: 



C= {0,1,...,A^-1}, axb=[a + b]^ 


for w = n: \F,fj:,^. ] = , . ,1 =1 


for w = n + 


rpabc 

[a+b+c]j^_ 

fpab 


= g»#"(^+c-[''+c]jv) 

[a+b]^[b+c]^ 

— g«#c([a-c]jv+6-[a+6-c]^) 

[a-c]jv[a+fc]jv 




Mab = e'—"' 




9, = e*2-f 



These anyon models describe some Chern-SimonsAVZW theories, e.g. SU(A^)i, 
for which the corresponding anyon models are Z^^ ^^^'^^ for A^ odd and Z^^^ 
for A^ even; and U(l)fc, for which the corresponding anyon models are Z^^/^^ for 2k 
even and for 2k odd. 

Of course, for Abelian anyon models such as these, each physical quasiparticle 
excitation has a specific anyonic charge and all fusion channels are uniquely de- 
termined, so superpositions of anyonic charge are not actually possible, but such 
models might occur as a subset of a non-Abelian anyon model, in which case su- 
perpositions of these charges could potentially occur. In any case, one may still 
perform interferometry experiments in these models to determine the charge of a 
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target anyon. Using b probes, we have: 

Pa = P2ao,b = 1^2^ + 2 \tirit2r2\ COS (^9 + ^j^ab^ + \rif \t2f (3.103) 
and 

PrA(K)= P(aJ-a-J){aJ-a-J) (3.104) 

Pt= Y: ^^"'^^^^^^^'i""'^^ \aj-a-f) {a'J-a';f\ (3.105) 

For zj,^^ with gcd(2w, iV) = 1 (i.e. the modular Z^r models), the charge classes are 
singletons = {a}, so a = a' in the fixed state density matrices. 



3.6.2 Fib 

The Fibonacci (Fib) anyon model (also known as SO(3)3, since it may be obtained 
from the 811(2)3 anyon model by restricting to integer spinsH] is known to be 
universal for topological quantum computation [|86il . It has two charges C = {l,e} 
and is described by (listing only the non-trivial F-symbols and i?-symbols, i.e. 
those not listed are equal to one if their vertices are permitted by fusion, and equal 
to zero if they are not permitted): 



C={l,e}, 1x1 = 1, 1 X e = e, exe = l + e 


[^fle/ = [Fself = 


" 0-1 0-1/2 
0-1/2 -0-1 


- ef 


J^ee g— i47r/5 J^ee ^i3n/5 


s= ^ 


1 
_0-l_ 




M = 


1 1 

1 -0-2 




di = 1, 4 = 0, V = ^(i) + 2 


9^ = 1, 9, = e^f 



where = is the Golden ratio. We denote the anyon model given by this with 
the complex conjugate values of the _R-symbols and topological spins as Fib. 

For b = e probes, we have Ci = {1}, C2 = {e} and 



Pi=Piii,e = I'tif 1^2!^ + 2 |tirit2T2| cos6' + \rif \t2f (3.106) 
P2=P7ei,e = ksl^ - 20~^ \tirit2r2\ COS + \rif \t2f (3.107) 

^ As a Chern-Simons or WZW theory, this is properly denoted as {G2)i, since SO (3)^ is 
only allowed for /c = mod 4. 
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Pta (1) = P(l,l;l)(l,l;l) + P{l,s;e){l,e;e) 



(3.108) 



p^' = P^{p(i,i;i)(i,i;i)|i,i;i)(i,i;i| 

+0"V(i,.;.)(i,e;.) |1, e; e) (1, e; e\} (3.109) 

Pr^ (2) = P(e,l;e)(s,l;s) + Pie,e;l)is,e;l) + P {e ,e;e){e ,e;e) (3.110) 
(P(e,£;l)(£,e;l) + P(e,£;e)(e,£;e)) 

X [\e,e; 1) (e,^; 1| + (3.111) 

We note that one can sometimes (approximately 69% of the time, when the target 
charge is not vacuum) perfectly distinguish the charges 1 and e with a single h = e 
probe measurement by setting the experimental parameters to: = \t2\^ = 1/2 
and 9 = n, which give pi = and p2 = 1 — ^ ^ .69. 

3.6.3 Ising 

The Ising anyon model is derived from the CFT that describes the Ising model at 
criticality |[8l. It has anyonic charges C = {1,0",^} (which respectively correspond 
to vacuum, spin, and Majorana fermions in the CFT). The anyon model is described 
by (listing only the non-trivial F-symbols and _R-symbols): 



C = {l,(j, V^}, lxa = a, a X a = 1 + a x ip = ip x ip = 1 


Fa Je 

r po-f/'o- 


(7 

aa L J (T(T 


\ 


1 1 

/2 V2 
1 -1 
/2 V2. 


ef 

= 

crw 

aa L ^ J 


= -1 


i?-- = e-'f, R^'' = e'^, Rf = R^'^ = e-'^ , Rf^ = -1 


S — ^ 


1 V2 1 
V2 -V2 
1 -V2 1 




M = 


1 1 1 
1 -1 

1 -1 1 






9, = 1, 9^ = e'^, 9^ = -l 



where e, / G {1, ip}. 
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For b — probes, we have Ci — {1, tp}, C2 — {cr}, and 



Pi — Pill,-!/' ~ Plp1pl,1p ~ Plt(>1p,1p ~ P-lpllp,'^ 

= |ti|^ Irsl^ + 2 ItirirstsI cos^ + |ri|^ [tal^ (3.112) 
P2=p;<,i,v, =P2ai>,i> = I^i|^k2|^ -2|iirir2i2|cos^+ Iril^l^al^ (3.113) 



PlA (1) =P(1,1;1)(1,1;1) + P{l,a;a){l,a;a) + P(1,V;^)(1,^;V) 

+P(V,1;V)(V',1;V') + P{'>l^,a;a){^,a;a) + (3.114) 



+P(V,V;i)(i,i;i) IV', "0; 1) (1, 1; 1| + P(V',V;i)(V,V';i) IV', "0; 1) (V', "0; i| 
~^\72 I-*-' ^' ^1 P(i,<t;<7)(^,<7;<7) |1, cr]a) {ijj, a; a\ 

+P(i,V;V)(i,^;^) ll,^':^') (1,V';V'I + P{i,V;V)(V',i;V') |1,V':V') (Va1;^/'| 

+P(v,i;v)(i,V';V') l'0, 1; V') (1, V'; V'l + P(v,i;v){V',i;V') IV', 1; V') (V', 1; V'l } (3.115) 

Pr^ (2) = p{a,l;<j){<T,l;a) + P(a,(T;l)(a,a;l) + P{<j,a;i>){a,<y;i>) + P{a,^;a}ia,^;a) (3.116) 



+P(o-,V';o-)(o-,i;a) \cr,'^;cr) (o", i;cr| + P(a,v;(T)(CT,v;a) |cr,V';'7) (f^, i;o-|) 

+P{a,a;f){a,a;f) k, 0"; V') (o-,0-;V'|} (3.117) 

For b — a probes, we have Ci — {!}, C2 — {cr}, C3 — {t/j}, and 

Pi=Pui,a = |ti|'k2|' + 2|tirir2t2|cos^ + Inl'ltsl' (3.118) 

P2 = p;.i,. = |ti|'|r2|'+|riHt2|' (3.119) 

P3 = Pi,^i,a = l^il^ 1^21^ - 2 |tirir2t2| cos6' + \rif \t2f (3.120) 

Pr^ (1) = P(l,l;l)(l,l;l) + P(1,<7;<t)(1,o-;<7) + P(1,V';V)(1,V';V') (3.121) 
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Pi = TsApr {p(i,i;i)(i,i;i) |i' 1; 1) 1; 1 



A ^ 

Pr^ 

1 , \ / 

H 7^P(1,<t;<t)(1,(7;(7) |1, CT', Ct) (1, CT] (T 



72^ 

+P(i,V;V){i,V;V') |1' V") (1, V'; V"!} (3.122) 



Pryl (2) =P((7,l;cr)(<T,l;(7) + P{cr,<T;l)(<j,cr;l) 

+P(o-,(t;?/))(o-,(t;?/)) + P((T,t/.;cr)(cr,i/.;cr) (3.123) 



^11 , N / , 

+ 2 (P('7,'^;1)(<7,<7;1) + P(<t,<t;V){<J,<t;V)) 

X [\a,a;l){a,a;l\ + \a,a;^){a,a;^\]^ (3.124) 

PfA (3) = P(i/>,l;i/))(V',l;5/') + Pi'ip,a-,a)(ip,a-,(T) + P(V',t/';l)(V',i/';l) (3.125) 

Ps = {pw,i;v)w,i;^) 1; ^) 1; ^1 



-p(^,^;i){V',V;i) IV^'V^; 1) (V^'^; 1|} (3.126) 



We note that one can always perfectly distinguish the charges 1 and ip with a single 
b = a probe measurement by setting the experimental parameters such that \ti\ = 
\t2f = 1/2 and 9 = n, which give pi = and ps = 1. 



4 Fractional Quantum Hall Double Point-Contact Interferometer 



After the detailed analysis of Section|3l one hopes that it has application in physical 
systems, and not just to the abstract idealizations that exist in our minds. In pursuing 
this hope, we turn our attention to fractional quantum Hall systems, since they 
represent the most likely candidates for possessing anyons and realizing braiding 
statistics of either Abelian or non-Abelian nature. 
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Fig. 3. A double point-contact interferometer for measuring braiding statistics in fractional 
quantum Hall systems. The hatched region contains an incompressible FQH liquid. Ss and 
Ds indicate the "sources" and "detectors" of edge currents. The front gates (F) are used 
to bring the opposite edge currents (indicated by arrows) close to each other to form two 
tunneling junctions. Applying voltage to the central gate creates an antidot in the middle 
and controls the number n of quasiholes contained there. An additional side gate (G) can 
be used to change the shape and the length of one of the paths in the interferometer. 

Indeed, a setup that is rather similar to the Mach-Zehnder interferometer described 
in Section [3] has been experimentally realized in a quantum Hall system [87]. This 
interferometer has, so far, only achieved functionality in the integer quantum Hall 
regime (though, even there, the physical observations are not completely under- 
stood 118818911 '). but it should be able, in principle, to detect the presence of braid- 
ing statistics [19019 11921 . and even discern whether a system possesses non-Abelian 
statistics {93]. Unfortunately however, there is a crucial and debilitating difference 
between the FQH Mach-Zehnder interferometer of [87J and the Mach-Zehnder in- 
terferometer described in Section [3l because of the chiral nature of FQH edge cur- 
rents, one of the detectors and its drain are unavoidably situated inside the cen- 
tral interferometry region. As a result, probe anyons accumulate in this region, 
effectively altering the target anyon's charge. This effect renders the interferom- 
eter incapable of measuring a target charge, and hence, useless for qubit readout in 
topological quantum computation. 

Fortunately, there is another type of interferometer that can be constructed in quan- 
tum Hall systems which z^' capable of measuring a target charge: the double point- 
contact interferometer. Moreover, such interferometers, which are of the Fabry- 
Perot type ||94ll , involving higher orders of interference, have already achieved 
experimental functionality in the fractional quantum Hall regime [|22ll . The dou- 
ble point-contact interferometer was first proposed for use in FQH systems in 
Ref. [95 1, where it was analyzed for the Abelian states. It was analyzed for the 
Moore-Read state [,31.[ . the most likely physical realization of non-Abelian statis- 
tics, expected to occur at z/ = 5/2 and 7/2 filling fractions, in Refs. [,96.69.97 .98,1 . 
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See also [I99I100I101I102I103I1041 ?1 for related matters. It was further analyzed 
for arbitrary anyon models, and specifically for all the Read-Rezayi states [36J, in 
particular, the one expected to occur at u = 12/5 filling fraction, in Ref. [85] (and 
subsequently analyzed for the = 12/5 Read-Rezayi state with homoplastic tech- 
niques in Refs. 11105110611 '). In all of these previous analyses for non-Abelian states, 
the results were given to lowest order in the tunneling amplitude, and only for target 
anyons that were assumed to be in a state of definite anyonic charge (i.e. already 
collapsed). In what follows, we provide expressions including all orders of tunnel- 
ing, both to explicitly display the unitarity of the quantum evolution and to account 
for potentially measureable corrections. Furthermore, we allow the target to be in 
a superposition of different anyonic charges, and relate the results to the analysis 
of Section [3l so that we now have a proper description of the measurement col- 
lapse behavior for these interferometers. Experimental efforts in realization of the 
double point-contact interferometer have been carried out for Abelian FQH states 
[|22I23I24I25I26I271| . Whether or not these experiments have conclusively demon- 
strated fractional statistics of excitations in the Abelian FQHE regime remains a 
topic of some debate I1107I108|F] . 



The double point-contact interferometer consists of a quantum Hall bar with two 
constrictions (point-contacts) and (at least) two antidots, Ai and A2, in between 
them, as depicted in Fig. [3l The constrictions are created by applying voltage to 
the front gates (F) on top of the Hall bar; by adjusting this voltage, one may 
control the tunneling amplitudes ti and t2. In the absence of inter-edge tunnel- 
ing, the gapped bulk of the FQH liquid gives rise to a quantized Hall conduc- 
tance: Gxy = I / {Vd^ — VsS) = ve^ /h, where the current through the Hall bar 
is J = {Id^ — Is^)- At the same time, the diagonal resistance vanishes: R^x = 
{Vd^ — VsS) // = 0. Tunneling current between the opposite edges leads to a 
deviation of G^y from its quantized value, or equivalently, to the appearance of 
Gxx oc Rxx 7^ 0. By measuring the diagonal conductance Gxx, one effectively 
measures the interference between the two tunneling paths around the antidot. The 
tunneling amplitudes ti and ^2 rnust be kept small, to ensure that the tunneling 
current is completely due to quasiholes rather than composite excitations. Treating 
tunneling as a perturbation, one can use renormalization group (RG) methods to 
compare various contributions to the overall current. Such analysis shows that in 
the weak tunneling regime, the tunneling current at a single point-contact has the 
dependence / oc where s is the scaling dimension/spin of the correspond- 

ing fields/anyons II28I101I102L It follows that the dominant contribution in this 
regime is from the field with lowest scaling dimension, which, in FQH systems, is 



^ One of the reasons for the uncertainty in interpreting the results of the experiments test- 
ing the Abelian statistics in the FQH regime is the fact that the statistical angle and the 
conventional Aharonov-Bohm phase acquired by a charged quasiparticle in a magnetic 
field are not easy to tell apart (this point is discussed in Refs. II109I95 1). From this perspec- 
tive, a non-Abelian FQH state might have an advantage, being that its effect from braiding 
statistics dramatically differs from the charge-background field contribution. 
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the fundamental quasihole. It should be noted that the quasihole tunneling is ac- 
tually relevant in the RG sense, which, in more physical terms, translates into the 
tendency of these point contacts to become effectively pinched off in the limit of 
zero temperature and zero bias. On a more mundane level, the quantum Hall liq- 
uid can be broken into separate puddles by the introduction of a constriction due 
to purely electrostatic effects (such as edges not being sufficiently sharp). In this 
regard, the recent experimental evidence [1 10], indicating that it is possible to con- 
struct a point-contact for which the u = h/2 state persists in the tunneling region, 
is reassuring. 



The two antidots are used to store two clusters of non-Abelian quasiparticles, Ai 
and A2 respectively, whose combined anyonic charge is being probed. The reason 
for two antidots, rather than just one (as has been previously suggested in [|95I96I97I98I85II ), 
is to allow for the combined target to maintain a coherent superposition of anyonic 
charges without decoherence from energetics that become important at short range. 
In particular, the energy splitting between the states of different anyonic charge on 
an antidot is expected to scale as (where L is the linear size of the dot) due 
to both kinetic (different angular momentum) and potential (different Coulomb en- 
ergy) effects [98 1 . On the other hand, for two separated antidots, this energy differ- 
ence should vanish exponentially with the distance between them, with suppression 
determined by the gap [|69ll . 



In order to appropriately examine the resulting interference patterns, we envision 
several experimentally variable parameters: (i) the central gate voltages allowing 
one to control the number of quasiholes on the antidots, (ii) the perpendicular mag- 
netic field, (iii) the back gate voltage controlling the uniform electron density, and 
(iv) a side gate (G) that can be used to modify the shape of the edge (and, hence, 
total area and background flux within) the central interferometry region. The rea- 
son for proposing all these different controls is to be able to separately vary the 
Abelian Aharonov-Bohm phase and the number of quasiholes on the antidots. In 
fact, having all these different controls may turn out to be redundant, but they may 
prove beneficial for experimental success. 



The target anyon A, is the combination of the anyons Ai, A2, and all others (in- 
cluding strays) situated inside the central interferometry region. In general, any 
edge excitation qualifies as a probe anyon, but since tunneling is dominated by the 
fundamental quasiholes, we can effectively allow the probes to have definite any- 
onic charge b equal to that of the fundamental quasihole. Letting (1, 0) and (0, 1) 
correspond to the top and bottom edge, respectively (also denoted as s =<-,->, re- 
spectively), the unitary evolution operator for a probe anyon B entering the system 
along the edge is given by 
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u 



rlr*e''^RABWABRcB 



i-lnrWBA 



(-1 + |r2|' Wab) Rcb nr^e^'^^BcRBAWBA 



(4.1) 



when the C anyons (those outside the central interferometry region that are entan- 
gled with A) are in the region to the right of central, where we have defined 



Wab = E (-tlt2e'^''^'''^RBARAB 



n=0 



[l + tlt^d^'^+'^^^RBARAB 



(4.2) 



The phases Oi and On are respectively picked up from traveling counter-clockwise 
along the top and bottom edge around the central interferometry region, and include 
the contribution from the enclosed background magnetic field. We note that when 
higher order terms are significant, it might be the case that tunneling contributions 
from excitations other than the fundamental quasiholes (which have different tun- 
neling amplitudes) are also important, but nevertheless proceed with considering 
all orders of tunneling in this manner. The tunneling matrices are 



(4.3) 



with j = 1, 2 for the left and right point contacts, respectively. We can perform a 
similar density matrix calculation as for the Mach-Zehnder interferometer, except 
with more complicated diagrams in this case. Sending a single probe particle in 
from the bottom edge (s =^) (which is effectively done by applying a bias volt- 
age across the edges), and detecting it coming out at the bottom or top edge gives 
the same form for the resulting density matrix as in Eq. (13.281) . except with more 
complicated ^ that are determined by using U of Eq. (14.11 ) for VU in Eqs. (13.91 - 
13.111) . To order \t\^ (for ~ \t2\ small), we find 



^ 1 - - \t2f - \ht2\ {e'^Mat + e-'^M:,,) (4.4) 

and 



+ \r,\' tit;e-*(^'+^")M:,, + Inl^ |t2|' 
^ \tif + \tit2\ {e'^Mab + e-'^Ml,,) + [tal' M^b (4.5) 
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where we have defined P = arg |t^t2e*'^^'^^"''|. We see that 

P7a'e,b + Paa'e,b - ^eb + ^2^ ■ (4-6) 

Here we have 1^2 1 as the probability of the probe B passing between anyons A and 
C, rather than , as in the case analyzed for the Mach-Zehnder interferometer, 
because of the location of C . The values for the two outcome probabilities (i.e. the 
e = 1 terms) to all orders are 



Paal,b — 51 ^ab 



J I |2 I |2 



izK^j^ ''';^^£,._„ . (4.7) 



dc 


In 


2 


1^2 


2 


dadb 1 + 


\ti\ 


2 


1*21 




1*1*21 


COS [P + 271 (Sc - Sa- Sb)] 



~ 1 - - \t2f - 2 |tit2| Re {e'^Mab} (4.8) 

Paal,b = 1 " Paal,b- (4-9) 

These are also the values of pla'e,b to orders when Meb = 1, but in general p^a'e.b 
does not have such a nice form. One may also perform edge current tunneling cal- 
culations employing the conformal field theoretic description of the edge modes in 
order to determine the effects of the source-drain voltage, the separation length be- 
tween the two point-contacts, and the temperature [95] (see also ^2811011102 103 106ll). 
The result of such considerations is essentially an interference suppressing Q-factor 
[recall the discussion following Eq. (13.201) 1 that decreases (with modulation) as any 
of these three quantities increase. There may be additional sources of interference 
suppression, such as switching noise [99| or edge-bulk tunneling H 1041 71. Though 
the suppression factor is only close to Q = 1 in certain regimes, we will again 
ignore it, but keep its existence in the back of our minds. 

As before, the target system collapses onto states with common values of p^ai 
generically producing a density matrix with non-zero elements that correspond to 
difference charges e with M^b = 1 (and Mab = Ma'b)- To first order, the behavior 
is essentially identical to that of the Mach-Zehnder interferometer which we pre- 
viously obtained, but the higher order terms may require more stringent conditions 
for superpositions to survive measurement collapse than just indistinguishability of 
monodromy scalar components (since this only guarantees proper matching to first 
order). Specifically, for superpositions of a and a' to survive, they must have 

^ jyc ^ / _ ^ ^ ( ^4 YO) 

c °'da\9aj c da' \9a' ) 

for all n, and some much more cumbersome condition for the survival of coherent 
superpositions corresponding to difference charge e. However, it seems that this 
condition is often equivalent to indistinguishability of monodromy scalar compo- 
nents for models of interest. In order to have p"^^^ b = 0, i.e. producing sometimes 
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perfect distinguishabilitjO], we require equal tunneling probabilities |ti| = \t2\, 
and cos [(3 + 27i (s^ - Sa - Sb)] = -1 for all iV^^ 7^ 0. In Eq. (II8OI) . we obtained 

an estimate for the total number of probes, N > (^j^^ , needed to collapse and 

distinguish a superposition of two anyonic charges in the target, with some level 
of confidence 1 — a. Translating this into the amount of time r necessary for such 

a measurement, we get the estimate r > (^T^^ ' where /tot is the total edge 
current. 

From these results, we find that when the target is in a state of definite charge a (or, 
more generally, in a fixed state with a E C^), the longitudinal conductance will 
be proportional to the probability of the probe injected along the bottom edge to be 
"detected" exiting along the top edge: 

Cxxcxp^,, , ~ |ti|V|t2|' + 2|M2|Re{e'^M,b} (4.11) 

which is exactly Eq. (7) in Ref. ['851. This is an experimentally measurable quantity, 
found by measuring the voltage between and D^. Using the side gate (G), one 
can vary P and, from the resulting modulation in the conductance, determine the 
amplitude of Mat- Indeed, the measurement of this quantity may be used to help 
properly identify the topological order of an unknown physical state. 

4. 1 Predictions for FQH States 

The results of this section are applicable to any FQH state. Because of their relative 
significance, we will only give the explicit details here for the Abelian hierarchy 
states, the Moore-Read state, and the (particle-hole conjugate of the) k = 3, M = 1 
Read-Rezayi state. The anyon models of these FQH states may be easily described 
in terms of those given in Section [X6l The application of this section's results to 
the entire Read-Rezayi series and the NASS states of [44J may be found in [[791 
(and partially in (851). hi all these FQH examples, it is important to remember that 
the electric charge and anyonic charge are coupled through Abelian terms for FQH 
states. Consequently, superselection of electric charge only permits superposition 
of anyonic charges that correspond to the same electric charge (i.e. are composed 
of the same number of fundamental quasiholes). 

4.1.1 The Abelian Hierarchy States [u = n/m) 

The Abelian fractional quantum Hall states can all be constructed from Z^r models. 
The general formulation in terms of K matrices may be found in UlllL but we will 

^ One can never have always perfect distinguishability for this interferometer, since it must 
be in the weak tunneling limit, which prevents ever having \tj\^ = 1/2. 
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describe the Laughlin and hierarchy states HI 1211 1311 1411 1511 that occur at filling 
fractions v = njm (with m odd and n < m). As shown in fST], the statistical factor 
of the fundamental quasihole in these states is 6* = ^ where p is odd and = 
1 mod m (which uniquely defines p modulo 2m). It follows that these states are 
described by 1^2m of Section |3.6.1[ in which a fundamental quasihole has anyonic 
charge [l]2m and electric charge e/m, while an electron (which has electric charge 
— e) has anyonic charge ["^Jg^. Using h = [l]2m probes, we have 



PaaO,b — 1 



ri ^2 



1 + iMale' 



Hif + \t2\ +2\tit2\cos[p + n 



2 Tip 



m 



(4.12) 



(4.13) 



and all permissible states are automatically fixed states with the target system's 
anyons A and C having definite charge 



Pa = l«.c; [a + cj^J {a,c; [a + c] 



2m I 



(4.14) 



4.1.2 The Moore-Read State {u = 5/2, 7/2) 

The anyon model corresponding to the Moore-Read state, expected to describe the 
u = 5/2, 7/2 plateaus, is given by 11791 : 



MR= Ising X Zl''"> ^ 

C = {(1, [2m]s) , (a, [2m + 1]^) , (^, [2m]s)} (4.15) 

(for m G Z), the restriction of the direct product of the Ising and Zg^^^^ anyon 
models to the charge spectrum C C Cisms x C„(i/2) in which the 1 and tp Ising 

charges are paired with the even sector of Zg^''^'* and the a Ising charge is paired 
with the odd sector of Zg^^^-*. Writing a = (aj, a^) G C, where aj G Cising and 
flz £ C (1/2), this is more explicitly given by 



C = {(1, [2m]g) , (a, [2m + l]g) , (V^, [2m]g)} , N^,, = Kl^N^lt. 


[^^ \ef- r'^' 


ei/i L '^^ J ez/z ' 


-ppab -pciibi -pCLibi 


pab _ paibi paz6z 
-"'c ^ci ^cz 


Sab — ^/^SaibiSa^bx 


Mab = Ma,b,Ma^b^ 


4 = (ia,, P = 4 





where all the symbols (A^, F, R, S, M, c/, 6) labeled by subscript 1 charges are those 
of the Ising model given in the table of Section I3.6.3[ and those with subscript Z 
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charges are those of the Zg^''^^ anyon model given in the table of Section [3.6.1[ The 
factor of \pl arises in front of the product of S-matrices because of the restriction of 
the charge spectrum. The fundamental quasihole has anyonic charge (a, [l]g) and 
electric charge e/4. The electron has anyonic charge (V^, [4]g). 

The target system's anyons A and C may by thought of as composed of n and m 
fundamental quasiholes, respectively (with electric charges nej^ and me/4), where 
n, m are integers (possibly negative). If n is even, the target's total anyonic charge 
may be in some superposition of (1, [n]g) and [rijg). If n is odd, the target anyon 
has total anyonic charge (a, [rijg). The same holds for the target's entangled partner 



C. We will employ the shorthand Xn = [x. 
that are fundamental quasiholes, 6 = (Xi, this gives 



[n]g), with X = l,a,tjj. For probe anyons 



Pi 

n 



P^ 

n 



Pa„a„lQ,<Ti 



I |2 I |2 



l + |t^t2|e<^+"?) 



\tif + \t2f + 2\tit2\cos (^P + rij 

I |2 I |2 
Fl ^2 



1 - Itital e 



I |2 I |2 /1 , I . . |2 



1 - (-1)^ 1^1^21^ e*2/3 



+ |t2|'-2|tit2| 



l + (-l)— COS (2/3) 



(4.16) 

(4.17) 
(4.18) 

(4.19) 
(4.20) 

(4.21) 



With these probes (and electric charge superselection), all charges are distinguish- 
able, so the charge classes are all singletons and interferometry will collapse any 
superposition of charge in the target onto a definite charge state. Of specific note is 
that for n odd, the interference is suppressed. The leading order modulation occurs 
at fourth order in t, and has the twice the modulation frequency. In fact, higher or- 
der harmonics enter as modulations in 2jP that are 4j*'^ order in t. When m is even, 
the probabilities and fixed states are 



P^a{^1„) P(ln,lm;ln + m){ln,lm;ln + m) ~l~ P(ln ,'/'m ; V'n + m) (In jl/"™ ; V'n + m) 



(4.22) 



Pl^l„ ~ T) ^ ^ 1 P(ln.lm;ln + m)(ln,lm;ln + m) l^™; ^m] l?l+m) (in; ^m] Ir 

1} (4.23) 
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P^A {/^an) — P{a„,lm;(rn+m){<r„,lm;Tn+m) P{(Tn,lf'm;Tn+m){Tn,i'm;(rn+m) (4.24) 



1} (4.25) 

)(V'n,lm; 

(4.26) 



1} (4.27) 

and when m is odd, they are 

Pr^ ('^In) = P(ln,0-m;0-„+m)(ln,0-m;0-n+m) (4.28) 

Pk;i„ "\72 '^"+"^) (in; '-'^mj '-''n+ml (4.29) 

Pr^ ('t(T„) P(o-n,0-m;lri+m)(o'n,0"m;ln + m) ~l~ P(o"n, jCTm iV'n + m) ("^n jO"™ i V'n + m) (4.30) 

(4.31) 

Pr^ ('^Vn) ~ P(V'n,iTm;crn+m)(V'n,CTm;cr„+m) (4.32) 

Pk^^ = IV'n, (T'm; (T'n+m) (V'n, CTjn] (7n+m\ (4.33) 

If one had a way to effectively suppress the tunneling of fundamental quasiholes, 
then the next most dominant contribution to tunneling comes from excitations with 
anyonic charge h = (1, [2]g), which are Abelian, and give (with different values of 



Px n 



= 1 



I |2 I |2 

ri r2 



2 (4.34) 



^ l^il^ + \t2f + 2 \tit2\ cos (^P + . (4.35) 

The Ising charges are obviously indistinguishable when the probe has Ising charge 
1, so superpositions of 1 and ip will not be affected by these probes (i.e. all al- 
lowed target density matrices are fixed states). If we had sufficiently good precision 
and control over the experimental variables to set them exactly to = |t2| — t 
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and cos + n|j = —1 for a = (1, [n]g) and b = {a, [l]g), then we would find 
Paai b = to all orders, providing a method of suppressing tunneling of fundamen- 



tal quasiholes. (These settings would givep^^^ 
this value of /3 for b = (1, [2]g) probes gives 



for a = {ip, [n]^).) Using 



Plnlnl0,l2 



1 - 



2n2 



1 - t2e™f 



:2t'^ [1 - cos (nn/A)] . 



(4.36) 
(4.37) 



Again, tunneling of these probes will be suppressed when n = mod 8, and 
the next most dominant tunneling contribution will be from b = [2]g), giving 



4.1.3 The Read-Rezayi State (z/ = 12/5) 

The anyon model corresponding to the state expected to describe the = 12/5 
plateau is [1791: 

RR3,i = Fib X ZSJ, (4.38) 

the particle-hole conjugate of the A; = 3, M = 1 Read-Rezayi statelZI- Writing 
) G C, where ap G Cpjj; and a-i G C„(3), this is more explicitly given by 



C = {(1, [n],,) , (e, [n],,)} , N^, = NZ^N^,^ 


\ef- l^* 


rpo-lfiici rpab rpapbp rpa^bz 




Sab SaphpSa-^bx 




da = dap, D = ^ 


'lO (0 + 2) 


da ^ap^az 



where n G Z, and all the symbols (A^, F, R, S, M, d, 6) labeled by subscript F 
charges are those of the Fib given as the complex conjugate of those in the ta- 
ble of Section U.6.2l and those with subscript Z charges are those of the Z^g^ anyon 
model given in the table of Section 13.6. 1[ The fundamental quasihole has anyonic 
charge (e, [1]iq) and electric charge e/5. The electron has anyonic charge (1, [5]^q). 
Being a direct product of Fib and an Abelian theory, universal topological quantum 
computation could be achieved through braiding quasiholes of this system. 



The target system's anyons A and C carrying total electric charges ne/5 and me/5, 
respectively, have az = [n]iQ and cz = [mlio, superpositions of 



° In general, the anyon models describing the k, M odd Read-Rezayi states may be written 
neatly as the direct product of anyon models RRfe,M = SO (3)^ x ^2(fcM+2)^^~^^''''^'' 
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ofjCf = l,e. We employ the shorthand (similar to before) Xn = {x, [n\-^Q), with 
X — 1,8. For probe anyons that are fundamental quasiholes, h — Si, this gives 



vX^e, - l^il' + + 2 |tit2| COS (/? - ny) (4.39) 



Pr„.„io,.i - l^il' + l^2|' - 20-2 \ht2\ cos (/3 - ny) . (4.40) 



With these probes (and electric charge superselection), all charges are distinguish- 
able, so the charge classes are all singletons and interferometry will collapse any 
superposition of charge in the target onto a definite charge state. We note that when 
the target has Fib charge e, the interference is suppressed, though still second order 
in t. The probabilities and fixed states are 



Pi ~ -r, I \ {P(ln,lm;ln+m)(ln,lm;ln+m) I"''"' ^ra'-, Iji+m) (Iri) Imj ln+m| 

"1"^ P(l„,£m;£n+m)(ln,£m;£n+m) 1-'-"' £n+m) (In) £n+m|} (4.42) 



(4.43) 



■''tsn ~ "P ^ ^ P(£n,lm;£n+m)(en,lm;£n+m) Imj £n+m) (^n; Imj £n+m| 

^ [l^n; ^mj In+m) (^ru ^mj ln+m| ~l~ l^n; ^n+m) (^n; ^mj ^n+ml]^ (4.44) 



By varying /3, one can distinguish whether ap, the Fib charge of a target anyon, is 
1 or £, without needing to know the precise value of the phase involved, because 
the interference fringe amplitude is suppressed by a factor of .38 for e. 

We emphasize that this provides the RR3 i state with a distinct advantage over the 
Moore-Read state with respect to being able to distinguish the non-Abelian anyonic 
charges that would be used in these systems as the computational basis states for 
topological qubits (i.e. 1 and ^ for MR vs. 1 and e for RR3 1). 
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